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ABSTRACT 

> 

Tlj" , Instead of imposing the Schrodinger equation to obtain the configuration 



, space propagator {x,t2\y,ti) for a quantum mechanical nonhnear sigma model, 

we directly evaluate the phase space propagator (a;|exp (^—^HAt^ \p) by ex- 
Q>^ I panding the exponent and pulling all operators p to the right and x to the left. 

I Contrary to the widespread belief that it is sufficient to keep only terms linear 

in At in the expansion if one is only interested in the final result through order 
Q I At, we find that all terms in the expansion must be retained. We solve the com- 

binatorical problem of summing the infinite series in closed form through order 
At. Our results straightforwardly generalize to higher orders in At. We then 
^ ' include fermions for which we use coherent states in phase space. For supersym- 

I metric A^=l and N=2 quantum mechanics, we find that if the super Van Vleck 

determinant replaces the original Van Vleck determinant the propagator factor- 
izes into a classical part, this super determinant and the extra scalar curvature 
term which was first found by DeWitt for the purely bosonic case by imposing 
the Schrodinger equation. Applying our results to anomalies in n-dimensional 
quantum field theories, we note that the operator ordering in the corresponding 
quantum mechanical Hamiltonians is fixed in these cases. We present a formula 
for the path integral action, which corresponds one to one to any given covari- 
ant or noncovariant H. We then evaluate these path integrals through two loop 
order, and reobtain the same propagators in all cases. 
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1. Introduction. 

In this article we settle a long standing problem, namely the direct computation of the phase 
space matrix element {x\ exp (^—^HAtj \p) for quantum mechanical systems with a kinetic term of 
the form a{x)p^ + b{x)p + c(f). The standard example, on which we concentrate, is the nonlinear 
sigma model describing a point particle coupled to an external gravitational field. Another example 
is the Hamiltonian for a point particle minimally coupled to an electromagnetic field. We shall also 
consider the inclusion of fermions, in particular super symmetric quantum mechanics. Although 
our results generalize straightforwardly to quantum field theory, we shall restrict our discussion to 
quantum mechanical systems. 

Through the work of Alvarez-Gaume and Witten on chiral and gravitational anomalies [1,2], 
and later the work of Bastianelli and van Nieuwenhuizen on trace (Weyl) anomalies [3,4], it is known 
that anomalies of an n-dimensional quantum field theory can in general be written as configuration 
space path integrals for nonlinear sigma models. However, the precise relation between the quantum 
mechanical Hamiltonian, which corresponds to the regulator in the quantum field theory and whose 
operator ordering is unambiguously fixed by this quantum field theory [5], and the corresponding 
action S'conf to be used in the configuration space integral, is unclear. For ehiral anomalies, the 
precise form of this action is not needed, as due to their topological nature the naive classical action 
5^,1 is sufficient. However, for trace anomalies the extra terms by which S^onf differs from S^-i are 
crucial. We therefore decided to settle the relation between S^onf "^ci in general. 

Before going on, let us eliminate a possible source of confusion. We are aware of an enormous 
amount of literature on configuration space path integrals, in which geometric operators like the 
Laplace-Beltrami operator and the Lichnerowicz operator play a role. In these articles these oper- 
ators are selected because of their nice mathematical properties. In the present article we do not 
begin by selecting such nice operators on mathematical grounds; rather, on physical grounds we 
begin with certain Hamiltonians H whose operator orderings are fixed a priori by corresponding 
quantum field theories, and our aim is to find the corresponding S'^.Q^f which produce exactly the 
same propagator as H. When it will so happen that these iS^onf do not have nice mathematical 
properties, so be it. For us, the relevant question is: which S^^^^ corresponds to which H? Of 
course, this question only makes sense if we also give a precise definition of how to evaluate the 
path integrals. We shall only consider perturbative expansions, and the precise rules are given in 
[1,3,4]. They will be reviewed in the appendix. (For the evaluation of anomalies, perturbation 
theory gives exact results). 

In the configuration space path integrals considered in [1,2,3,4], the internal symmetry sector 
or the fermionic sector were still treated by an operator formalism. We shall treat the bosonic 
and fermionic sectors on equal footing, starting with an operator formalism and deducing from it 
the corresponding path integral formalism. Our final path integrals will be genuine path integrals, 
without sectors where one still uses operators. To achieve this, we need to evaluate the matrix 
element (x| exp (^—^HAt^ \p) for bosons or {x,fj\ex.p (^—^HAt^ \P:0 fermions where {fj\ and 
\0 are fermionic coherent states. 
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At first sight the task seems appalUng (which may account for the absence of this result in 
the literature) because expanding the exponent one finds momentum operators p in all possible 
places between position operators x which appear in the arbitrary metric g.ij{x). Pulling all p to 
the right and keeping track of all possible commutators of p and x seems a hopelessly complicated 
combinatorical problem. It seems universally believed that it is sufficient to expand exp ^— ^.ffAt^ 
into 1 — ^HAt, evaluate the matrix element {x\H\p) = h{x,p){x\p) or {z*\H\z) = h{z* , z){z*\z) 
and reexponentiate to obtain exp (^—^h{x,p)At) {x\p) or exp (^—^h{z*, z)At) {z*\z) (here 1^) denote 
bosonic coherent states) [2,5]. This we have found to be incorrect: one must keep all terms in the 
expansion of the exponent, but for each term one needs to take into account only a finite number 
of commutators between p and x if one is only interested in the result to a given order in At. This 
makes the problem tractable after all. Including fermions, one obtains in addition to x and p also 
operators and V'^. Again, expanding the exponent, all terms contibute but only a finite number 
of anticommutators of and need be kept in each term, and the transition element (also called 
propagator) (x,f/| exp (^—jrHAt^ \PiO is obtained to any given order in At. 

To avoid confusion, let us define what we mean by 'to a given order in At'. We say that the 
propagator is given up to order / in At if it can be written as a product of some function f(x,y; At) 
and a factor [1 + ©(At)'"*"^/^]. The function / itself will in general be singular in At. 

The (incorrect) neglecting of terms of order (At)^ in nonlinear sigma models should not be 
confused with the (correct) neglecting of terms of order (At)^ in the approach of Feynman and 
others for Hamiltonians of the form T{p) + V{x), where the Trotter formula[6] actually proves the 
correctness of this procedure. Nonlinear sigma models are simply not in the class of models to 
which the Trotter formula can be applied. 

It is often claimed that after integrating out the momenta, one can use the resulting propa- 
gators {x.j\cxp (^—j-HAt^ to order At to obtain the configuration space path integrals by 
taking products of these building blocks and integrating over the intermediate points x^ . . . 
However, there is a problem: usually one wants to use local actions S^onf, and the matrix ele- 
ments {Xj \ exp ^—^HAt^ l^j-i) cannot in general be written as the exponent of a local action. 
For example, for the nonlinear sigma model we consider below, one finds terms AtR{x) and 
R^j{x){x — yy{x — yy, and whereas the former can be written to order At as / dtR{x), the latter 
term cannot be written in such a way. In fact, the nonlocal {xj\ex.p (^—j-HAt^ l-^i-i) satisfy the 
product property of path integrals, and thus trivially produce in a path integral the propagator 
(x| exp (^—j-HT^ \y). The problem is rather to find a local action S'^onf' which also produces in a 
path integral this propagator. From the example of the point particle in curved space it is known 
that S^onf contains local terms of higher order in h, but in general the problem of constructing S^^^^ 
from a given Hamiltonian (with given operator ordering) seems unsolved. DeWitt has recently made 
great progress in this direction by establishing a formal (unregularized) connection between H and 
'S'conf using the concept of time-ordering [7]. Using his results and the two and three-loop results of 
[4] on trace anomalies, we shall give the exact local action S^oni which corresponds to any, covari- 
ant or noncovariant, Hamiltonian H with at most two p operators. We have not been able to give 
a rigorous proof of this result, but we have checked our result through two-loop order, and give 
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general arguments which hold at higher loop order. 

Actually, Feynman already suspected that for nonlinear sigma models integrating out the 
momenta would not lead to path integrals with the classical action in the exponent [8]. In the 
literature an indirect method has been used to obtain the propagator {x,t2\yTti) in configu- 
ration space. Namely, one requires that it satisfies the Schrodinger equation H^{x,t2\y,ti) = 
ih (5/9^2) (3J)*2ly)*i)) where is the Hamiltonian in the ^-representation, {x\H\x') = H^6{x—x'). 
In this way, DeWitt found that for small t2 — , this propagator does not approach Feynman's 
result, but an extra term proportional to the scalar curvature is present [9]. We shall be able to 
trace its origin by following the direct calculation step by step. Bastianelli and one of us[4] found 
the extension of DeWitt 's result to fermions. They took a very general ansatz for the action S^^^f 
in the configuration space path integral, and after evaluating this configuration space path integral 
by a loop expansion, they were able to determine the propagator by again imposing the Schrodinger 
equation. 

The reader may start wondering at this moment why we use at all the direct but cumbersome 
method of evaluating J dp {x\ exp ^—^H^ \p){p\y), rather than just following the heat kernel ap- 
proach. We have two reasons. First of all, since Feynman, most physicists derive path integrals as 
we do, and it is desirable to understand this approach in detail. Secondly, it will help us in deter- 
mining the action S^^^f. In fact, in flat spacetime, one obtains directly S^onf after integrating out 
the momenta which are introduced by (x| exp (^—^^H^ \p); only in curved spacetime complications 
set in. 

There are two distinct results one is interested in: the propagator and its trace. Let us 
first consider the case of flat spacetime. Then both results are obtained from a formula which is 
obtained in the spirit of Feynman by inserting complete sets of x and p states or coherent states, but 
in the former case the propagator can only be written as a phase space path integral provided one 
adds boundary terms depending on a complex classical trajectory, while for the trace one directly 
obtains the path integral without the extra boundary terms and without the complex classical 
trajectory. These boundary terms can be found in textbooks[5], whereas the results for the trace 
formula were used by Alvarez-Gaume in [2]. In section 3 we explain the difference between both 
approaches. In curved spacetime, as wc already mentioned, even for small At the propagator may 
contain terms which cannot be written as local actions, and in these cases the straightforward 
Hamiltonian approach does not yield the action ^conf • As we already stated, we have a formula for 
the correct S^^^^f in general, and it contains the same kind of boundary conditions as found from 
the Hamiltonian approach in flat spacetime. 

Since it is sometimes claimed that phase space path integrals are less well-defined than config- 
uration space path integrals, and that path integrals for fermions cannot be as rigorously defined 
as for bosons, and because we were also ourselves confused for a long time, we decided to give a de- 
tailed account of the resolution of these problems, starting from scratch in section 2. In particular, 
we want to stress here that path integrals for fermions are phase space path integrals, and perfectly 
well defined in terms of coherent states. (The Grassmann integration makes them actually more 
convergent than their bosonic counter parts). This is very well explained in some textbooks[5]. 
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and we have nothing new to add, but for completeness we also include these aspects in section 5. 
In reference [5] the Bargmann-Wigner formalism of holomorphic functions is used to define inner 
products, Hilbert spaces, and only afterwards one tackles quantum mechanics. This we found an 
unnecessary detour; one can directly define coherent states in quantum mechanics, and use the in- 
ner product for these coherent states which follows straightforwardly from the (anti)commutation 
relations of the annihilation and creation operators. 

The fermionic path integrals need Grassmann numbers ry" and 77" (a = 1 . . . n) which are 
independent. However, the Hamiltonian which corresponds to the regulator Ipip of spin-^ quantum 
field theory, contains only one kind of fermionic operators, (One identifies the Dirac matrices 
7" with \/2'0"(f)). This satisfies the equal-time anticommutation relations = (5"^, but 

in order to obtain creation and annihilation operators, one needs operators and V"" satisfying 
l^a^ ^afj _ ^ah ^^^^ {t/j"- , tp^} = tp''^} = 0. It is Suggested in the literature that one may double 
the number of fermion species in the Hamiltonian, and replace the spin connection oj,ig^iiipfipi by 
with a = 1, 2. Afterwards one is supposed to divide the final results by 2"/^ in order to 
undo this doubling. This practise is well-known for one- loop problems in quantum field theory [10], 
but it is incorrect to use it in this context because the one-loop anomalies of a quantum field 
theory correspond to multiloop corrections of the corresponding quantum mechanical model[3,4]. 
We resolve this problem by adding a second set of free fermions. 

We begin in section 2 with a discussion of phase space path integrals based on x and p eigen- 
statcs. In section 3 we discuss the bosonic phase space path integrals based on coherent states. 
Section 4 contains the heart of this paper: a direct evaluation of the propagator through order At of 
bosonic nonlinear sigma models. We use here x and p eigenstates. The result leads us to the com- 
plex classical trajectories, and the trace formula. In section 5 we study the N=2 supersymmetric 
nonlinear sigma model and discuss the relation between operator ordering of the Hamiltonian and 
the supersymmetry charge on the one hand, and the various supersymmetric or non-super symmetric 
actions on the other hand. In section 6 we repeat the analysis of section 4 including fermions and 
using fermionic coherent states. In both cases, the resulting propagator can be written as a product 
of three factors: the classical action, the (super) Van Vleck determinant, and a term which involves 
only the scalar curvature and which is related to the trace anomaly in two dimensions. In section 
7 we discuss the N=l supersymmetric nonlinear sigma model, in particular the resolution of the 
problem of fermion doubling alluded to above. Again the propagator can be written as the product 
of the super Van Vleck determinant, the exponent of the classical action, and a term involving 
the scalar curvature which determines the trace anomaly of in this case a spin-^ field. Section 8 
contains the conclusions. In particular, we give here, without proof, the one to one correspondence 
between H and S^^^ for any, covariant or noncovariant, H. In the appendix we perform the 2-loop 
calculations based on configuration space path integrals, using covariant and noncovariant actions 
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2. Phase space path integrals in general. 



Path integrals were (almost) introduced into quantum mechanis by Dirac in 1933 [11]. In 
those days, quantum mechanics was obtained from classical mechanics by replacing the canonical 
brackets of the Hamiltonian formalism by corresponding quantum commutators. Dirac set out 
to find a formulation which was based on the Lagrangian in order to be able to retain manifest 
Lorentz invariance at the quantum level. His starting point were the following relations of classical 
mechanics between coordinates, conjugate momenta and action 

^S,,{x,y-t)=p{t) ■ ^S,,{x,y;t) = -p{0) (2.1) 

where S^i{x,y;t) is the classical action for a point particle as a function of its initial coordinate y 

and its final coordinate x. These relations follow immediately from the Euler-Lagrange variational 
t 

equations for jdt' L{x{t'),x{t')) at fixed t. Dirac found a corresponding operator equation in 



quantum mechanics as follows. Introducing "moving frames" \x,t) = exp ^^Ht^ \x), where \x) 
are the eigenstates of the position operator x = x{t=0) at time t = 0, and \x, t) are the eigenstates 
of the position operator x(t) = exp (^^Ht^ xexp (^—^Ht^ at time t, he noted that the inner product 
between these complete sets of states, called "transformation function" by him, satisfies a relation 
similar to (2.1). Namely, defining a function U{x,y;t) by 

{x, t\y, t=0) = exp (^U{x, y; tU (2.2) 



/d 
dx' {x, t\p(t)\x', t) (x', ils, (=0) = -m-^ {x, t\y, (=0) 

X (2-3) 



one has 



—U{x,y;t) j {x,t\y,t=0) 

where p{t) is the momentum operator at time t. A similar relation holds for p{0). Dirac then 
introduced an operator U{x{t),x{0)) whose matrix elements are the overlap functions in (2.2) 

{x, t\U{x{t), xiO)) \y, t=0) = U{x, y; t) (2.4) 

He called such operators "well-ordered", meaning that all x{t) should be put on the left and all x(0) 
should be put on the right, where they can be replaced by their eigenvalues x and y, respectively. 



^These moving frames project a ket \ip) onto the corresponding Schrodinger wave functon ipix, i) = {x, t\il}), 
and in this x representation, p{t) is represented on tpix, f) as = —ih-^, as follows from {x, t\p{t)\x' , t) 

■ — x'). The overlap function in (2.2) gives then the probability amplitude for the point particle 

to be at position x at time t, if it is initially at time t = at position y. In the nonlinear sigma 



models we consider we shall use as inner product \J g{x)i^x\x') = 5{x — x') where we define 5{x — x') 
by Jdx' f{x')5{x — x') = f{x). There is then an ambiguity in the representation of p{t), namely p^ = 
—ihg^-^g^" with g = detg^j{x), which we fix by requiring that p{t) be hermitian with respect to this 
inner product. The result is a = —1/4. 
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Combining (2.3) and (2.4), Dirac obtained operator equations corresponding to (2.1) 

He then noted that according to the completeness postulate of quantum mechanics one has 

{x, t\y, t=0) = / dxi . . . dxjv_i {x, t\xj^_i,t]\f_i){x]\f_i,tj^_i\x]sf_2, tN-2) ■ ■ ■ ^il^j t=0) (2.6) 



and stated that for fixed tj and small h only stationary points of U contribute. 

The relation (2.5) suggests that the operator U is the kind of quantum mechanical extension 
of the classical action Dirac was looking for. In fact, he went further and carefully stated that the 
propagator {Xj_^i,tj^-^\xj,tj) "corresponds to" the classical action S^i{Xj^i,Xj;tj^-^ — tj) for small 
tj j^i — tj. He did not say that they are equal but suggested that they were proportional, in which 
case he obtained an equivalence principle for small h, provided the stationary trajectories of U 
coincide with the classical trajectories of S^.^. 

In this article we consider quantum mechanical systems with a kinetic term of the form a{x)p^ + 
b{x)p + c{x), where x and p have n components (x* and with i = 1 . . . n), also called nonlinear 
sigma models. The models we shall consider are point particles, or point spinors coupled to an 
external gravitational field. An example of the kind of Hamiltonians we study is given by 

H = \g-'/\x)n,g'/\x)g^^{x)rr^g-'/\x) - ^K6cd('^(^))(V'aV'^) iWp) (2-7) 

where g{x) = det g^j{x), a and /? run from 1 to 2, and a,b and i,j run from 1 to n. The Riemann 
tensor is defined in (5.4), while 

^i=Pi + ^.^iabi£)rJ^ ; {raA'p} = ns'^%p ; [Pi,x^=-in5i (2.8) 

This is the Hamiltonian for N=2 supcrsymmctric quantum mechanics. For the A^=l case tpf = 
■02 = "^^/) can be identified with the covariant derivative + \(^iabl"'l''' if o^^c identifies 

with Vh^"". In fact, for the A'^=l case, and also for the purely bosonic case, the particular 
ordering of all operators in H (including the position operator x which we shall write from now 
on for notational simplicity without hat), is dictated by the n-dimensional quantum field theory 
from which they come. We now explain this point; however, readers not interested in quantum 
field theory may skip the following paragraph and simply restrict their attention to the particular 
order ings of H which we consider. 

There is a general construction [12] for any quantum field theory which yields the regulator 
TZ one must use in order to obtain consistent anomalies. If the Jacobian for an infinitesimal 
transformation of a classical symmetry is given by 1 + J, the anomaly is given by 

A= lim TrJexp(^) (2.9) 



6 



For a scalar field in n dimensions this construction yields TZg = g~^/^dj^g^/^g'^Wjg~^/^ for the con- 
sistent regulator which preserves general coordinate invariance. On the other hand, the consistent 
regulator which preserves Weyl (local scale) invariance but, as a consequence, breaks Einstein 
(general coordinate) invariance, is given by d^g^^'^g'^Wj. We call these regulators consistent be- 
cause they lead to anomalies which satisfy the consistency relations. The corresponding quantum 
mechanical Hamiltonians are obtained by replacing by ^p^. The matrix elements of in the 
^-representation will not be needed (they are given in the footnote above), since all we will use 
are the commutation relations and the fact that Pi\p) = Pi\p) on eigenfunctions \p). However, for 
completeness we mention that in the ^-representation with inner product (/i,/2) = / dxy/gf*f2 
one must replace by jQ^^^^d^g^^^ . Then the quantum mechanical operator corresponding to 
the regulator becomes g~^^'^dig^^'^g^Wj, which is indeed hermitian with the inner product with y/g. 
For spin-^ fields TZf,,.,. = g^^'^Ip Ip g^^^'^ which can be written as 'Rf^r = ^/er + i-^ where Tl^fer 
the first term on the right hand side of (2.7) and R is the scalar curvature. The Dirac matrices 
7" satisfy {7^,7^} = 2r/"'' and are represented in the corresponding quantum mechanical model by 
operators V2i^j satisfying {ipj,i^^} = ff^- In this case there is thus no ordering ambiguity in the 
Hamiltonian. For other applications [2] one may need an TV = 2 action with (1, 1) supersymmetry 
in which case a and /3 range from 1 to 2. In this case we shall fix the operator order ings in another 
way, because there is no quantum field theory for which the regulator corresponds to the N=2 
Hamiltonian. The result is that the last term in the Hamiltonian in (2.7) is then the same (up to a 
sign) as the 4-fermion term in the classical action. Even though there are no ordering ambiguities in 
the Hamiltonian which corresponds to TZg or TZ , it remains to be seen whether the corresponding 
action we will produce in the configuration space path integrals is the invariant classical action. 
We shall obtain the anomalies of the corresponding n-dimcnsional quantum field theory from the 
propagator {x,t\y,t=0) by multiplying with the Jacobian J(x), putting x = y, integrating over 
X and taking the limit t tending to zero, retaining only the finite part, (t^^ plays the role of the 
square of the regulator mass, M^). As a result the classical terms (of order h~^) in the propagator 
as well as the one-loop (^-independent) terms will not contribute, and only the two-loop terms of 
order h will contribute for small ^ in n = 2 dimensions. This is not in disagreement with Dirac's 
results, since he only made a statement about the propagator for h ^ at fixed t. 

Dirac did not write down an equation relating the propagator {x,t\y,t=0) for small t to the 
classical action S^i{x,y;t) because he probably did not know the normalization factor between 
them. Before him, it was already known to some mathematicians (for example Wiener in 1923[13]) 
that one could define functionals of functions x(t) by time discretization and taking ratios of two 
such functionals. For example, for the functional = Jdn exp (^J^V{x)dt^ with V{x) a 

local smooth function of x{t) such as x'^{t), and Jd/i = jDx exp (^—^m Jq x'^{t)dt^ , the following 
limit exists and defines the symbol / Dx 



F[V[x]] 
F[V = 0] 



lim 

AT-* 00 



F[V,N] 
F[0,N] 



(2.10) 



Here F [V, iV] contains the time discretization 

f I 2 

F\y, N] = dxi . . . dx^ exp I ^ --■m[xj_^i - Xj) + V{xj)e 



(2.11) 



with Ne = t and Xq = x{t=0) = y is fixed. The path integral in the denominator can be viewed 
as a measure for the path integration. In the Hamiltonian approach, to which we now turn, one 
derives the measure from first principles. This has the advantage that one need not worry whether 
the measure in (2.10) is the only reasonable choice of measure. 

Path integrals as such, and their normalization, were found by Feynman in 1948[8], who went 
back to Dirac's work, but inserted not only the A'^ — 1 complete sets of intermediate coordinate 
eigenstates, but also N complete sets of momentum eigenstates. He considered quantum mechanical 
systems with a Hamiltonian of the form H = T{p) + V{x), and using the relation 

exp (^-^HA?j = exp ^-^TAtj exp (^-^VA?j + 0{Atf (2.12) 

he obtained phase space path integrals for a point particle with n components by substituting the 
following expression into (2.6) 

= ^(x^+i|exp (^^TAt^ exp (^^VAt^ \xj) 



Pj+i 

J (27r^)" 



exp (^i^j+i - Xj)Pj+i^ exp (^-^T{pj+^)A?j exp (^-^Vixj)A?j 

(2.13) 

where At = tj^-^ — tj. For T{p) = i^p^ , he could integrate over p and obtain the configuration 
space path integrals of Dirac, but with their normalization 

(777 \ n/2 / ^ \ o 

^-^) exp (^-5(x,.+i,x,.;i,+i - t,)j + 0{AtY (2.14) 

Here S = T —V , but in Euclidean space one may start with exp {—y^HAt) in (2.12) and then one 
finds exp {-\Se) in (2.14) with Se = T + V. 

In the actual loop by loop calculations we perform in the appendix, we do not use (2.6) and 
(2.14), but rather go over to so-called mode- variables in the continuum theory a/,. One can carefully 
follow the steps leading from the discrete-time variables to the variables , and then take the limit 
N CO. The result is that one obtains a nontrivial measure d/i = (27re?i)^"' ^ ( fjH" ) d'^k- 
Using this measure, propagators and vertices follow, and loops can be calculated. For a detailed 
derivation of this measure see [3,4,14]. 

Since Feynman, path integrals have taken an enormous fiight. There are two approaches: the 
phase space approach and the configuration space approach. In the phase space approach there 
always comes a moment where exp (^^HAt^ \pj) has to be evaluated. It is our contention that 
the evaluation of this basic matrix element is always done incorrectly for nonlinear sigma models, 



8 



even though for models in fiat space the usual procedure gives the correct result. Namely, without 
exception (to our knowledge), it is always assumed that the following is a correct approximation if 
one is interested in a result which is corrrect up to order At 



1 exp (^-^HA?j \p^) = {xj\l - ^HAt\pj) (2.15) 



We have evaluated this matrix element exactly for nonlinear sigma models, and found, as conjec- 
tured by Feynman, that one cannot neglect terms of order {At)'^ in the expansion of the exponent 
for any k. Most of these At will be re-exponentiated and appear as —-^^p^At in the exponent; 
they yield the normalization factor in (2.14) after integrating over p. However, some terms remain, 
proportional to the same exponent times a polynomial in p. One should really keep all terms in the 
expansion of the exponent, move all p to the right and all x to the left, and each time a p passes 
a x it produces a factor —ifi. In terms of rescaled variables q = p-\fejh. where e = At one then 
obtains expressions of the form 

(x|exp (-^^At) \y) = {2'Kh)--g-^l\x)g-'l\y) {^^^'^ j d^q exp (^-^^q.{x^ - y^)^ 

exp {^^1i9^^i^)'i?j 1 + V^A(g,a;) + enB{q,x) + ... 

(2.16) 

where A contains one or three factors of q and contains the terms due to one p,x commutator, while 
B contains zero, two, four or six q momenta and contains the terms due to two p,x commutators. 
Integration over q replaces q by mix — y)/\f€E, and since x — y is of order -y/e, we can stop the 
expansion in (2.16) at B for a two-dimensional quantum field theory. The factors g~^^'^{x) and 
g~^^^{y) come from the inner product 

{x\p) = {2nn)-^'^g-^l\x) exp (^^^x') (2.17) 

which is a direct consequence of the orthogonality relations \/ g{x){x\x') = 5{x — x') and {p\p') = 
S{p — p'). Putting A = B = {) amounts to the approximation in (2.15), but is is clear from the 
way we have displayed the further terms in (2.16) that after integration over the terms with A 
and B will contribute corrections to the propagator of order \/lE and eh. This last term will give 
the trace anomaly for n = 2. For a quantum field theory in dimensional space (or spacetime), 
one must determine the expression within square brackets in (2.16) to order (e?i)""/^, which is 
cumbersome. In the configuration space approach[l,2,3,4] which uses the Schrodinger equation, 
one must compute "graphs" on the worldline with Uq loops, which is also cumbersome. Only for 
chiral anomalies where the Grassmann integration over the fermionic zero modes brings in factors 
of eh one need only expand the configuration space path integral to second order in quantum 
fluctuations [1,2]. There are then already enough factors of h to ensure that it makes no difference 
whether one takes the classical action or another action which corresponds to a different operator 
ordering of the Hamiltonian. 
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3. Path integrals with coherent states. 

In this section we want to discuss phase space path integrals based on coherent states. These 
are the path integrals one uses for fermions (because for fermions ^ and are more like a and 
then like p and x). There are some misconceptions about such path integrals, having to do 
with the observation that a path connecting a point in phase space to another point in phase 
space cannot in general be a classical path. We begin by quoting Schulman[15] who in his very 
readable textbook on path integrals already anticipated the extra terms in (2.16): ". . . No one (to my 
knowledge) has made a serious investigation of the neglected terms e{zj — Zj_i) [he refers here to the 
approximation in (2.15) and the extra terms in (2.16)]. My own guess is that they can contribute 
and that this contribution will be related to the operator ordering problems in quantum mechanics 
. . ." . We agree with his observation on operator ordering: if the exponent of the Hamiltonian of 
the nonlinear sigma model would have been well-ordered in Dirac's sense, no extra terms would 
have appeared in (2.16), but physics (the n-dimensional quantum field theory) has given us another 
ordering. However, Schulman considers at this point coherent states \z,t) (see below), rather than 
the eigenstates |x, t) and |p, t) we considered so far, and claims that the extra terms are produced 
because the curves z{t) are more singular than the Brownian motion paths that enter the usual 
path integrals (for which (Aa;)^ ~ At). We do not agree with him on this point; already in (2.16) 
which does not use coherent states one encounters the extra terms. And he goes on[15]:". . . A well 
known feature (see Section 31) of the phase space path integral is that paths in phase space are 
discontinuous so that a term eAjo need not go to zero any faster than e . . .". Indeed, if p{t) is the 
time derivative of x{t), and the paths x{t) are continuous but not differentiable, then the paths 
p{t) will be discontinuous. However, looking up this Section 31, one discovers that he considers the 
following path integrals 

(3.1) 

He then states[15]: ". . . we interpret the formula (3.1) as a sum over paths [correct]. The integration 
variables suggest that one is adding trajectories in phase space [incorrect in our opinion] . . . However, 
. . . the fact that a single point in phase space determines the classical path, causes some difficulty 
[indeed]. Here are two interpretations . . .". And he goes on to consider either nonclassical paths in 
phase space for which Xj is continuous but Pj jumps at the end of each time segment, or classical 
paths for which p{t ■ + ^At) = p- and x{tj + At) = Xj^^ with tj^i — tj = At. We claim that in 
(3.1) one is considering ordinary paths in configuration space, with beginpoint Xq = y and endpoint 

= X. For coherent states things are more subtle as we now discuss. 

We define coherent states for bosonic point particles by 

|2) = e«^^|0) ; (z*| = (0|e^*« ; [a,a^] = l (3.2) 

where z* = (z)* is the complex conjugate of z, and z = {x+ip)/\/2h while z* = {x—ip)/V2h. (This 
corresponds to a harmonic oscillator with H = \p^ + \SP' = Ti{a^a+\). For H = ^p^ + \muP'x'^, 
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one must use Xy/muj and p/y/muj in these definitions). The vacuum is annihilated by a, and one 
has the following inner product and decomposition of unity 

r dzdz* 

{w*\z)=e'^' ■ 1 = j^^-\z)e-' '{z*\ (3.3) 

The symbol ^^^^ always denotes ^5^, and we stress that z and z* arc not independent complex 
variables. The propagator {w* ,t=T\z,t=^) = (it;*|exp (^—-^HT^ \z) is then equal to 

{w*,t=T\z,t=0) = / rr -^—L{w*\expi--HAt) e""^-!"^-! 

^ 7=1 ^""^ \ ^ J (3 4) 

exp (^-^HA?j |zjv_2) • • • e-^i*^i exp (^-^HA?j \z) 

with At = T/N. Note that this is an exact result for any A'". In section 4 we shall compute for small 
At the matrix elements of exp (^—^HAt^. Anticipating this calculation, let us define a function 

h{z*^^,zj) by 

exp (^-^HA?j \zj) = exp z^)At j (3.5) 

Then the integrand in the propagator can be written as 

where z^ = z and z"^ = w*. This is the formula referred to in the introduction, from which the 
propagator and its trace are obtained. 

For the trace anomaly one should take the trace of the propagator and evaluate 

Tiexp(^-^HA?j = j'^^^e-'"'^w\At\z,t=Q) (3.7) 

where w* must be equated to {z)* . The integration in (3.7) is then again an ordinary Gaussian 
integral. To prove (3.7), one may either use (3.3) or expand the coherent states into harmonic oscil- 
lator eigenfunctions \n) = (n!)~^/^(a^)"'|0). The result reads then indeed X]„(n| exp i^—^HAt) \n). 
The extra term e"'^*^ in (3.7) completes the exponent of (3.6) to an action provided we identify 
z = Zq with ZjY, hence z = Zq = z^ 



exp / dt 



-z*z — ■^h{z* ,z) 



(3.8) 



This path integral has then periodic boundary conditions. In terms of x and p the kinetic term 
reduces to the familiar — ^px. For the computation of other anomalies one must evaluate the trace 
of (u;* I J exp (^—^HAtj \z), where J is the quantum mechanical operator which corresponds to the 
Jacobian of an infinitesimal symmetry transformation of the quantum field theory. Depending on 
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the problem, the presence of J can be translated into a change in boundary conditions (periodic, 
antiperiodic, etc. [2]). 

For other problems than anomalies, one may need the propagator itself. In order to still be 
able to write (3.6) as a path integral, wc no longer have the extra exponent due to the trace 
operation but rather we must add and substract a suitable term by hand. In principle, any term 
will do, but for practical calculations it is very useful to decompose z{t) and z*{t) into a classical 
part which satisfies the equation of motion, and a quantum part over which one still integrates 



with the measure ^ / e ^ ^. Terms linear in the quantum field then cancel from the action, and 



dz*dZi 

3 

2ni 

the corrections begin with the harmonic (quadratic) approximation, and then higher terms. To 
achieve this one considers two independent classical trajectories Zf.i{t) and iu*i(t), uniquely specified 
by Zf.i{0) = z and w*^{T) = w*. Of course, the value of Zf.i{T) at i = T is then not w, nor is w*^{t=0) 
equal to {z)* . Rather, to avoid confusion, wc keep denoting them by z^i{T) and u'*j(0). We then 
complete (3.6) by adding and substracting a term w*z^i{T). This yields as path integral 



1 

-z^z — —h{z*^z) 



(3.9) 



One can also treat the beginpoint and endpoint more symmetrically, by adding one-half the above 
term, and another term given by ^w*i{0)z. In this case one finds the path integral[5] 

gi[«,*.,,(r)+<,(o).] J ^ _1 _ -^^^ _ ^h{z*,z) ^ (3.10) 

The boundary conditions are now: z*{T) = w* , z*{<d) = ^^^(O), z{T) = z^i{T), z{0) = z. In the 
action one is then to replace z{t) by %(*) + z^^{t) and similarly z*{t) by iy*i(t) + ^qu(*)- The 
quantum fields vanish at the boundaries 

= V(^) = <u(0) = <u(r) = (3.11) 

and we still have z = {x + ip)/ \/2^, z* = {x — ip) j \/2^, with real x and p. The measure ^2tH^ ~ 
^2t^ defines an ordinary Gaussian integral. One could define a classical p and x by 

<i{t) = {x^i{t)-iPci{t))l^ ; ^ciW = (^ciW + ^PciW)/^ (3-12) 

In that case, x^i{t) = (w*l(t)+z^l(t))^/h/2 and in particular x^i{t=0), x^i{t=T) are complex (because 
?w*j(t=0) z*). Similarly for There is no problem with classical paths connecting two points 

in phase space because we have here two classical paths which do not connect w* and z, but one 
of which starts at z and ends wherever the dynamics may take it, while the other starts at w* and 
ends somewhere else. No discontinuity or jumps of paths need be considered. 

To evaluate (3.10) for the harmonic oscillator is trivial [5]. Terms linear in quantum fields 
cancel whereas the terms quadratic in quantum fields are independent of the classical fields and 
yield an overall constant, so that the classical factor in front of the integral in (3.10) contains 
the whole result up to a constant. This constant can be fixed by considering the limit T — >^ 0, 
or by doing the path integral over z and z* explicitly. For more complicated systems, one can 
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evaluate the path integral by expanding the action in terms of classical and quantum fields, but 
then z^^{t) is not an independent complex variable, but rather expressed in terms of (2:qu(t))* by 
Zquit) = - + (^qu(*))*- We can define x^^ and p^^ by z^^ = {x^^ + ip^^) /V2h and 

-^qu = (^qu ~ Wqu) /V^, but then also x^^{t) and Pqu(*) are complex. Using 

x^^{t) = x{t) + ^{z,i{t) + wl,{t)) ; p^^{t)=p{t) + ^^iz,,{t)-w;,{t)) (3.13) 

it is clear that the Jacobian for the transformation of the integration variables x{t) and p{t) to 
Xqu(i) and Pq^{t) is unity, so that the measure becomes 

Dx^^Dp^^ ; x^^{tj) e[-oo + aj,+oo + aj] 
Pquitj) ^[-oo + bj, +00 + bj] 

where Qj = ^ (^cil^y) +^ci(*j)) ~ h {^cSS-j) ~'"^ci(*j)) complex constants for each ty 

In other words, x^^{t) and p^^{t) lie on a line in the complex plane parallel to the real axis. In 
a loop expansion, the integrand is analytic in all variables, and we may simplify the integrand to 
^qu(*) = {xqu{t) + Wanit)) and z\^{t) = {x^^{t) - ip^M) l^f^ with real x^^{t) and p^^it), 

and the integration measure becomes 

Dx^J^p^^ ; -00 < ajqu < 00, -00 < < 00 (3.15) 

For fermions we consider operators and i^) satisfying {V-',V'^} = 1? and independent Grass- 
mann variables f] and r]. This will force us to double the number of fcrmionic operators, from il)^ 
to as in (2.7). The operators i/)" and V^"^" are the defined by + ii/^f )/y2 and (V^" - i%j:2)/\f^, 
respectively. We continue in this section with the case n = 1, omitting the superscript a. The 
coherent states are now defined by 

= e^^'^IO) ; (^1 = (016*?^ (3.16) 

satisfying V'l??) = ^|??) and = (rylry. The inner product and decomposition of unity read 

formally the same as in the bosonic case 

(^10 =e^^ ; 1 = Jdfjd^\Oe-'^^{fi\ (3.17) 

but now the ordering of df] and d^^ as well as the ordering of f] and ^ in f]^ is important, while the 
factor (27ri)~^ is absent from the measure because we now have Grassmann integration according 
to which Jd^^ = Jdfif] = 1. The results in (3.17) follow directly by expanding the exponent; one 
finds 1 = |0)(0| + |1)(1| with '0^|O) = |1), which is evidently correct. We can repeat the same steps 
as before, leading to (3.6). The propagator is then given by 




I 



(3.18) 
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or, in the symmetric case, by 



1 ■ . i 



(3.19) 



Again we obtain the well known extra boundary terms in these path integrals with fermionic 
coherent states[5]. On the other hand, the trace formula is given by 



TVi = (0|i|0) + (l|i|l) = j d^df?e^«(r?|i|0 



(3.20) 



where now ^ and rj are independent variables. The proof follows directly from Grassmann integra- 
tion and shows that stands to the left of df). The trace of the propagator becomes then 

N-1 / ^ N-1 \ 

j d^dfj Yl dfjjd^j exp I 77^ + fj^N^i - %_iC7v-i • • • - ViCi + ^i^-^Yl ^(^j+i' (^-^l) 

This can be interpreted as a phase space path integral with J dt {—fj^— ^h), provided one interprets 
^{t=0) = = —^{t=T). By putting the term fj$, in the exponent to the far right, one again finds a 
path integral, now with J dt {fj^ — ^h), provided 'fj{t=0) = —fj = —fj{t=T). Taking the symmetric 
case, we recover (3.8), but now with antipcriodic boundary conditions for ^ and f]. For other 
anomalies with J not equal to unity, the boundary conditons on ^ and fj may be different. For 
example, for chiral anomalies one obtains periodic boundary conditions for ^ and fj [2]. 



4. Transition amplitude for the bosonic case. 

We will consider the matrix element 

(x|exp(-^i?) |p) (4.1) 

with 

H = l9-'/%9'/'9'^Pj9-'/' - l^n'R (4.2) 

We are working from now on in Euclidean space; to obtain the results for Minkowski spacetime 
one should replace e by ie. For ^ = j^Ef) H is the regulator for a classically conformal invariant 
scalar quantum field theory in n dimensions [3]. Expanding the exponent in (4.1), we define 

{x\[h) |p)^^Af(x)p'(xb) (4.3) 

1=0 

where Af[x) is a c-number function and denotes a homogeneous polynomial of order I in the 
momenta. 
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In order to compute the transition amplitude {x\ exp (^—^H^ \y) to order e compared to the 
leading terms, it will turn out that we only need terms on the right hand side of (4.3) with I = 2k, 
2k — 1, and 2k — 2. We find, defining = g^-' {x)pj^Pj, 



2k 



1 



(4.4) 



Since this is the term containing the maximal number of p's, it can be easily computed because all 
p operators are just replaced by the corresponding c-numbers when acting on \p). 
The next term is 



2k-l 



—ihk- 



1 /I 



^P 



2 V 2 



— ih 



k-1 



i^^9'')PJ 



^p 



k-2 



(4.5) 



P]PkPl 



In this expression one of the ^'s acts as a derivative, whereas the other 2k — 1 are replaced by the 
corresponding c-numbers. The first term in (4.5) comes about when the derivative acts within the 
same factor H in which it appears, and is multiplied by k since there are k factors of H. The 
second term arises if this derivative acts on a different factor of H. For this to occur there are (2) 
possible combinations, and taking into account that there are two p's in each factor of H wc get 
an extra factor 2. Notice that in both cases the terms involving a derivative acting on g cancel. 
The last term we have to calculate is 



A2k-2ix)P 



2k-2 



K'k [ -p 



k-l 



^g'^ {d, logg) {dj logg) + ^g'^ {d,d^ logg) + J {d,g'^) {dj logy) 



^P 



32 

k-2 



2^ 



k-3 



^g^^(d,d,g'^^) + l{d,g^^) (a,/') 

+ \ {d,g^') {d,g^') + Id'" {d,d,g^^) 



Pjpl 



+ 



\g^' (d^g^'^) (d,g^^) + \g'^ (d^g^^) (d^g^^) 



PkPlPmPn 



2^ 



^5*^5-- (5,/') (d^gP'i) 



PjPkPlPmPpPq - o^fi k[ -p 



k-l 



R 

(4.6) 

In this expression two of the p's act as derivatives, except for the last term which is already of order 



p'2h 2 

The first set of terms appears when both derivatives act within the same factor H; again there 
are k terms of this kind. 

The next set of terms arises when only two of the factors H play a role. There are four 
possibilities: {i) one p from the left factor acts on the right factor, while another p from the right 
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factor acts within the right factor, (ii) the first p acts within the first H, while the second p acts 
within the second H, {iii) both p's come from the left H, but one of them acts inside the left H 
while the other acts on the right H, and {iv) both p's from the left H act on the right H. In all 
cases it is easy to see that again the derivatives on g cancel. 

The following set of terms comes from combinations using three factors H, hence its overall 
factor (3) . There are again four cases: (i) a p from the first H and a p from the second H hit the 
third H, (ii) one p acts inside the factor H in which it appears whereas a p from another H hits 
the remaining H (there are 3 terms of this kind), {iii) a p from the first H hits the second H, and 
a p from the second H hits the third H, and (iv) of the two p's from the first H one acts on the 
second, and one on the third H. 

Finally, the term with (^) involves four factors H, such that one p from one H hits another 
H, and the other p from one of the remaining factors H hits the last H. 

The reason further terms do not contribute can be most easily seen if we rescale q = \/\p- 
Then the transition amplitude becomes 



(x| exp \y) = J d^P {x\ exp \v)^\y) 



g-^'\x)g-^l\y){27:n)-^ (- ] J 



q exp I 



iix-yy 



eh 



(4.7) 



-1/2 



k=0 1=0 



where we have used (2.17). So only the ^2fc-i ^2k-2 terms contribute through order e compared 
to the leading term ylgfc- 

The sum over k in (4.7) can be performed, leading to 

{x\ exp( - ^h) \y) = g-''\x)g-y\y) {^^ne)-''" J cTq exp (-Ig'^q^Qj + ^^^^^^ 
h iV7h!^\ (d,g^^) qj - \g^^ [d.g'^') q^^q, } 

^9'' {d, log g) [dj log g) - ^g'^ {d,d^ log g) - ^ {d,g'^) [dj log g) 



+eh< 



(4.8) 



+ 



^,9V{l 



+ 
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Ikmmin 



QjQkQlQmQpQq 



We can now complete the square in the exponent and integrate out the momenta g^, since the 
integral becomes just a sum of Gaussian integrals which can easily be evaluated. The problem is 
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now to factorize the result such that it is manifestly a scalar both in x and in y (a 'bi-scalar') 
under general coordinate transformations. We expect, of course, to find at least the classical action 
integrated along a geodesic (see (4.11)). In the expansion of this functional around x{0) = x, we 
recognize many of the terms in (4.8). However, there are terms left over. They combine into R 
or R^j, while expansion of g{y) yields terms with dlogg or derivatives thereof. With this in mind, 
we write the result in a factorized form, where in one factor we put all terms which possibly can 
come from expanding some power of g{y), while into another factor we put the expanded action 
and curvature terms. It is quite nontrivial, and an excellent check on the results obtained so far, 
that this is at all possible. The resulting expression is 



(a;|exp( - -H)\y) = {l^Ue)-^/^ g-"\x)g-"\y) 



g^'\x) + g^/\x){y - xf [d,g^'\x)) + \g^'\x){y - xfiy - x^ (d,d^g^l\x)) 



exp ( -2^9ijix)iy - x)\y - xy 



1 fl 1 



1 1 



+ 2 (^4^ i^k9ij{x)) {y - x)\y - xy{y - x) 

~T2^ - l9mnix)r7;{x)rUx)j {y - x)\y - xy{y - xf{y - x)' 

+ {^i^i-Y^^^R{^)-Y^M^)iy-^ny-^y + o{e^i^) 



(4.9) 



where the Ricci tensor is defined by R^j = Ri^kj and Rij^i is given in (5.4). The terms within the 
first pair of square brackets are, through order e, equal to g^^^{x)g^^^(y) and cancel the factors 
g~^^'^{x)g~^^^{y) in front of the whole expression. Note that from the term in the exponent in (4.9) 
it follows that the difference {y — x) is of order y^, thus one indeed finds this expansion through 
order e. Similarly the terms within the second pair of square brackets are of order e or less. The 
terms with S^.g'j^- and its square are the first two terms in the expansion of an exponent. This 
suggests to exponentiate all terms, yielding 



{x\exp(^--Hj \y) = (27r?le) 



-n/2 



exp 



\9ijix) + ^d,^gij{x){y - x)'' 



+ ^ (dkdigijix) - ^g^^ix)r^ix)TUx) ) {y - xny - x) 



{y - xY (y - x) 



+^^{l^-^)Ri^)-^Riji^)iy-^yiy-^y + ^(^'/') } 



(4.10) 



All terms in the exponent except the last two just correspond to an expansion around x of the 
classical action, which is equal to the integral along the geodesic joining x and y of the invariant 
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line element (cf. [9]) 



t=-e 



dt^9iMt))x\t)x\t) 



\9ij{x){y - xfiy - xY + ^d^gijix)iy - x)\y - xY (y - xf 



+ 



{x)mx)TUx) {y - xYiy - xY {y - xf{y - x) 



+ -0{y 
e 



(4.11) 



when boundary conditions x{—e) = y, x(0) = x are imposed. 

Our final result for the transition amplitude can thus be written as 

1 



{x\ exp (-|^) \y) = (27r;ie)-"/2 exp 



2^ 12 



1 



-{R..(x) + Ri^iy)) {x-yfix 



{R{x) + R{y)) 

yY + oie^') 



(4.12) 



This shows that (4.10) is, to the order we are expanding, symmetric under the exchange of x and 
y. One may check that it satisfies the composition rule 

J dz,/^{x\exp (-|^) \z){z\exp \y) = (x|exp (-f ^) b) (4-13) 

A quick way to check this is to use normal coordinates, in which case one only needs to retain 
the leading term in the classical action (since dj^gji^ = d^j^djQ/^i^ = 0), while gij{z) = gij{x) + 
^Rj^^ji{x){z — x)''{z — xY through order e. Taking the opposite point of view, we can impose (4.13) 
as a consistency condition on the amplitude. This fixes the higher order terms in the expansion of 

(4.12) , with exception of the terms proportional to the scalar curvature, which are related to the 
trace anomaly. One should expect this to be the case since we can always add such a term to the 
Lagrangian and interpret it as an external potential, which should not be fixed by the requirement 

(4.13) . 

The transition element can of course also be evaluated using path integral methods. One should 
be careful however in choosing the correct action in the configuration space path integral, since 
our Hamiltonian (4.2) is not Weyl-ordered[16] or time-ordered [7]. The corresponding Weyl-ordered 
or time-ordered Hamiltonian is (see e.g. [17,7]), dropping a noncovariant term proportional to the 
product of two Christoffel symbols for reasons given in the conclusions. 



(4.14) 



Since it is -ffweyi that corresponds to the classical action in the path integral, the (Euclidean) action 
that corresponds to H in (4.2) is = ^ci ~^ l^'^ I -^j agreement with the results in [3]. 

We have obtained the result for {x\ exp (^—^H^ \y) from a direct computation; the extra R and 

Rj^j terms are due to ordering the p operators in all terms of the expansion of exp (^—^h\ The 
terms with R^j can be expressed in terms of the classical action as follows. By using 



D'/\x,y;€) 



-l'g'l\x)g^l\y) 



xfiy 



(4.15) 
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where D{x,y;e) = det D^^(a;, y; e) is the Van Vleck determinant [18,19] 

d d 

Dij{x,y;e) = y; e) (4.16) 

we can rewrite (4.12) as 

{x\ exp \y) = (27rn)--/2^V2 exp (-l^.i) [l + ^eh - ^) {R{x) + i?(y)) j (4.17) 

where D = g~^l'^{x)D{x, y; e)g~^^'^{y) and we have neglected terms of higher order in e. Note that 
D as well as S^i{x,y;e) transform as biscalars under general coordinate transformations {D itself 
is a bi-density [20]), hence also the infinitesimal transition amplitude (xjexp (^—^H^ \y) does not 
depend on the coordinates chosen. Obviously this result is not affected by the terms proportional 
to the scalar curvature. 

The Van Vleck determinant gives the harmonic approximation of the quantum terms in the 
path integrals [15]. For example, the Feynman factor (^)"^^ in (2.14) is a Van Vleck determinant. 
Since -^S^\{x,y]e) = —pj with Pj the momentum conjugate to , D is the Jacobian for the 
transition of the phase space measure dx2dpidxi to the configuration space measure dx2dydxi 
where y divides the interval (xj^, X2) into two smaller intervals of size TjAf. To the order we working 
in, D for the half-interval is equivalent to over the whole interval and reproduces the terms. 
One might think that one could explain the R^j terms in the propagator by taking as exp(S'(,oiif) 
the expression ^^/^ exp(S'ci) . However, this is incorrect because the path integral for this S^^^{ 
produces its own R^j terms, on top of those contained in D, as was already noted in [3]. 

We conclude that the final result can be written as the product of three factors: one term 
which is only related to the trace anomaly, the exponent of the classical action, and the square root 
of the Van Vleck determinant. 

As one would expect, the transition amplitude is the Green function associated with the 
Schrodinger equation (diffusion equation for Euclidean time) 



exp (-iij) Is) = ng-^'^(x)S{t)i"{x - v) (4.18) 
where 

H{y) = -\n^g-"\y)^^g"\y)g^={y)^ - ^e^^i?(y) (4.19) 
To prove (4.18) we begin with the identity 

(x| exp(-^F)y) = Idz ^/^{x\H\z){z\ exp{-^H)y) = Hix){x\ exp{-^H)y) (4.20) 

If we would expand the propagator around y, and act on it with H{x) as in (4.20), we obtain the 
same result as acting with H{y) on our propagator (which is expanded around x), because the right 
hand side of (4.18) is symmetric in x and y. Using the footnote in section 2 and (4.18) we may 
check (4.10). 

We can now simply compute the anomaly in n = 2 by taking the trace in (4.12) and singling 
out the e independent part. The result is[3,4] 

^^calar ^ -^R (4.21) 
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5. The classical and quantum Hamiltonians for the N=2 supersymmetric nonlinear 
sigma model. 

As an example of a ease where the operator orderings in the Hamiltonian are unambiguously 
fixed by requiring that the corresponding configuration space path integral has the original classical 
action in its exponent, consider N=2 supersymmetric quantum mechanics. The classical action in 
Minkowski time is given by 



This action can either be obtained by dimensional reduction of the A^=(l, 1) supersymmetric non- 
linear a model from 1 + 1 dimensions to one time dimension[l], or directly by requiring invariance 
under (1, 1) supersymmetry transformations with parameter e^. One finds either way the classical 
supersymmetry transformations, and the classical supersymmetry charge. The latter is given by 



To define the corresponding quantum generator, it is helpful to introduce fcrmions with tangent- 
space indices, because their brackets are field- independent. Therefore define ■i/'^ = e"(0)Xa where 
^iA^ab = 9ij- Then the action in Minkowski time becomes 



(5.1) 



and is invariant under the following supersymmetry transformations 



= i (eiX2 - £2X1) = e^cr2X' 
,5x' = -i0V2e-rV'/''x*^ 



(5.2) 



(5.3) 



Rijabi^) = diUJjab + ^ia^'^jcb - (i ^ j) = Rijki(r)e^ei 



(5.4) 



The supersymmetry charge becomes 




(5.5) 



So far, this is standard supergravity. 

To quantize this system, we first observe that there are second class constraints 



(5.6) 
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The Poisson brackets are given by {p{'ipa)''4''^} p = "'^"^'^a/Ji where p(V'S) = (*^® minus sign 

in the Poisson bracket is not a matter of convention but follows from requiring the Heisenberg field 
equations to hold). Using [C^,C^^p = —i5^^5^p^ we find the Dirac brackets 

1 , ^^■'^'> 

X0.0 X 

which is as usual half the Poisson brackets. Thus the quantum anticommutator reads (we set ^ = 1 
in this section) 

|C,V'j|=5"'5a/J (5.8) 



Furthermore, the conjugate momentum of (f)^ is given by 

P^ = 9^34^ + \^^abrai^^ ^ [ft, = -^^l (5-9) 

We now choose as quantum supersymmetry charge 

Q^^-^ = eiicl>)ry/\g-y' ; ^. ^ p. - la;,„,C^^ (5.10) 

We choose this particular ordering because then gi'^''^ is hermitian, and hence also i^^'^") will be 
hermitian. To prove that this qL""^ is hermitian, one may show that after hermitian conjugation one 
gets three extra terms which cancel due to the "vielbein postulate" that the total covariant deriva- 
tive of the vielbein (with Cristoffel and spin connection) cancels. We evaluate the quantum anticom- 
mutator {Qa''\Qf''^} and show that the result is equal to S^^H^'^''^ with by definition the 
quantum Hamiltonian. Classically, the Hamiltonian is given by ^g^-''^i'^j — ^Rabcdi^){'^a''Pa){'^'^'^'^)-> 
but quantum mechanically there are many possibilities for the Hamiltonian, because {tp^^ V'^l^'^of 
vanishes classically for any X, but equals S'^^S^^X'^^ at the quantum level. We define the Hamil- 
tonian as the square of the supersymmetry charge because, as we shall discuss in section 6, when 
this is used in (a;,?7|cxp (^ — eH^ ly^x) it gives the same propagator as the path integral in 
which the original supersymmetric action appears. Hence, wc have chosen what might be called 
super symmetric ordering for Q^a^"^ and H'^'^^\ Also, with this ordering presciption, Qq^"'* and H^'^"^ 
transform as scalars under general coordinate transformations applied to operators (cf. [7]). 
The direct evaluation of {Qa'^\Q^^^^} is tedious, but straightforward. Using 

one finds, keeping all operators in the order they come, 

{qI^''\ Q^r^} = g'^' {Kp^'''^>A'^, - V'^C K, vr,-] - i + a;,-,'vr,] g-^l' (5.12) 

(We used the vielbein postulate 9^6^ = — r*^e^ — oj^^^e^ and found that all F terms cancel). The last 
term vanishes classically, but at the quantum level it is nonvanishing; in fact, it just covariantizes 
the TT^ term to the complete d'Alembertian for a scalar particle! Using 

W'^^=-\^iiabrj>l (5.13) 
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and the cyclic identity 



''abed 

''abed 

we see that also the four-fermion term is proportional to d^p. Our final result reads 



(5.14) 



(5.15) 



Thus, imposing the quantum algebra {Q^^''\ Q^^"^^ } = ^a^H^"^"^) with given Q^a^^ has fixed the 
operator ordering in as given in (5.15). Note that this quantum Hamiltonian is already Weyl- 

ordered the term with tt^ yields upon Weyl-ordering again a term see (4.14), whereas the 

Rijj^ term contributes — cancelling the bosonic contribution. To show this latter result, write 
the Weyl-ordered operator corresponding to |-Ra&c<iV'a'0aV'|a'0^ = 5-Ra&c(iV'"'0^V''^V''^ (using (6.3)) as 



o ^abed 



eit„ t,)d{t„ Qe{te, uWi^'r^" - e{t,, t,)9{t„ tMtd, t,)ri^'^''r + 22 more (5.16) 



where t„ denotes the time associated with ip"', and the limit tg^,ti„t^,t^ — t is understood. We 
do not have to include the curvature term in the Weyl ordering because it commutes with all the 
fermions. To compute the difference between (5.15) and (5.16), we decompose unity in a sum of 
products of ^-functions, and write the operator in (5.15) as follows 

^R,,,ari^'r^^ = ^it;„,»V'^V''' [o{ta, h)0ih, Qe{t,, t^) + e{t,, t,)9{t„ t^)e{t^, t,) + 22 more " 

(5.17) 

Now substract (5.17) from (5.16), using the anticommutation relations (5.8). When we make use 
of identities such as 

+ e{t^, tj9{t^, t^)em, ta) + ^(ie, ta)e{t^, t^td, h) + oit^ h)9{h, tjeit^, t^) 

(5.18) 

then, in the limit that all times coincide, we find that the difference equals —^R- So the two 
extra terms due to Weyl ordering the two terms in in (5.15) cancel, and indeed the corre- 

sponding configuration space path integral which follows from this particular ^(i*^) has the classical 
supersymmetric action in its exponent. 



■tWe obtain this covariant result if we use a Weyl ordering in which we consider tt as an independent 
operator, and not p. Note that DeWitt in [7] at this point switches notation, and denotes g^^x^ by Pj (his 
equation 6.7.26). Weyl ordering (which he discusses when p and tt coincide) in terms of his new p and x 
then indeed agrees with our observation that (5.15) is already Weyl-ordered. However, from a canocical 
point of view, tt is not an independent operator. We discuss this point further in the conclusions. 
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6. Transition amplitude for N=2 supersymmetric sigma model. 



We will now compute the matrix element 



(x,?7|exp (6.1) 



where 



H = \9-"%g'"9'^^,9-"' - \n'Ra,ai {r.i'i) {mi) 

T^i=Pi- Y^iab'^a'^i ; a = 1,2 

As shown in (5.15), this is the quantum Hamiltonian for the N=2 supersymmetric nonlinear sigma 
model. The states |x) and (^| are eigenstates of the operators and ■0'' respectively, V'^lx) = x"lx) 
and {fjl'tjj^ = {fj\ff'. They are coherent states, \x) = exp (■;/;" x") 1 0) and {fj\ = (0| exp(r7"'0") , and 
satisfy the completeness relation 1 = J(^d^d^^ . . . d^'^d^,'^) exp(— ^"^") |^)(^|, and inner product 
{v\0 = cxp(f?^). 

Compared to the appearing in the regulator for the trace anomaly of a spin-^ field {Ttfer = 
we have twice as many fermionic degrees of freedom (a = 1,2). The case with half 
as many fermionic degrees of freedom (with a = 1 in (2.8)), does correspond to a quantum field 
theory for spin-^ fields and will be discussed in section 7. We define 

Analogously to the bosonic case we expand the exponent in (6.1) and define 

2k 

{x, r?| {H)''\p, O^Yl ^'(^' ^' ^) ^\P' (6.4) 

1=0 

Again we will only need the terms with Z = 2A;, 2A; — 1, and 2k — 2. It will be convenient to express 
the Hamiltonian as = a + /3 + 7, with 



P = -ing'iu,,,r^'g'l%g-^'' (6.5) 



I 2 "jao/T r 2 

The leading term is the same as in the bosonic case 

k 



The next term equals 

fc-i 



Blk-^{x,f),Op"'-' = ^2fe-i(^y'-' - ikh 9''u;,,bf]''epj (6.7) 
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The first term comes from k factors a, and is identical to what the purely bosonic case yields, 
whereas the second term arises when we combine k — 1 factors a with one factor and replace all 
operators by their corresponding c-number or Grassmann number values. The last term we have 
to compute is 



2k-2 



k-2 



+ Ikh' Qp^j g^^{d,logg)u;^,,te - lk{k - l)n' l^^p'j g'^{d,g'')cv^,,tePkPi 

- Ikh' (^P')' ' [9-'/'d,{g''9'^'^jab)v''e + {9''u;,^u;^^^ 

- ^k{k - i)n' ' ' 9''u^iab9'''^kcd{v^('^'"' - n^n^e^^jPi 

(6.8) 

There are four possible combinations that give a contribution: [i) k factors a, yielding the same 
expression as in the bosonic case, (ii) k—1 factors a and one factor (3; now one p, either from one 
of the a's or from /3 acts as a derivative, again either on one of the a's or on /3, yielding the next 
four terms, (in) k — 1 factors a and one factor 7 and replacing all p's by their c-number values, and 
(iv) k — 2 factors a and two factors (3, which gives rise to the last term. Notice that to evaluate the 
last two terms we first needed to normal order the fermionic operators which produced extra terms 
proportional to S"'^ according to (6.3). In detail, we write 'ijj'^tlj^'tp'^il^'^ as '4)^ip'^5^'^ — 'ijj""tp'^tlj^ip'''. We 
can now proceed in the same way as in the bosonic case, by doing the sum over k, extracting a factor 
exp(— 2^p^), and performing the Gaussian integral over the momenta. In the last term of (6.8), 
the integral over pjPi contributes a term with gji which cancels the complete lo'^ term in the line 
above it, plus a term with (y — x)j{y — x)i. In this result, the contribution proportional to uj'^ff^'^ 
is the square of the second term on the right hand side of (6.7) and is removed by exponentiating. 
The two-fermion part of the last term in (6.8) survives and is also exponentiated. We then compute 
the transition amplitude 

(x,f?|exp (-^h) \y,x) = J tTpcT^'cT^ e-«"«(x, f?| exp (-|^) b, x) (6-9) 

Using the inner product 

{x,fj\p,0 = {x\p){m = (27rn)-/V'/'(^)e*^"e*^« (6.10) 
we find for the amplitude 

(a;,r?|exp (-^i?) \y,x) = (27re?i)-"/2 exp (^-^S^ - ^5^) 



1 _ ^ehRix) - ^Rij{x){y - x)\y - xY + \enR^{x)n'' 



(6.11) 
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where 



Sb = -^9ij{x){y - xfiy - xy + ^4(?i (a;)(y - xfiy - {y - xf 



+ ^ - \g^^{x)TT^{x)TUx)) {y - xfiy - xy{y - x)\y - 

is the expansion through order e of the length of the geodesic joining x and y (cf. (4.11)), and 

~Sf = -nS^btx' - Hy - xyu;,,,{x)tx' - \n{y - xy{y - xy (d,u^M + oo,,'^{x)oo^,,{x)) ftx^ 

1 



(6.13) 

If we add an extra boundary term ^'5aftV'ci('-')V'ci(0) = ah^'^ '^\\^) (6.13), then the sum of the 
first term in (6.13) together with this extra term becomes the leading term in the fermionic action 



~Sp = Sp- U^^ff^lm ■ Sp = h^dt (^d^^ri;" + x'u,,t,ri;' - ^hR^.^r^'ri^'^'^ (6.14) 

— e 

We can easily check that the expansion through order e of (6.14) indeed equals the expression in 
(6.13) when the equations of motion arc imposed. These read (of course we also need the bosonic 
part of the action to find the full equations of motion) 

x' + r^.x^x'^ - R'ja.x^r^' + \ng'^ {djR,,^) r^'r^^ = 

^« + - W = 

We can now expand the Lagrangian in a Taylor series around its value at t = 0, and then do the 
trivial time integrations. This yields 

S = eL(0) - ^e'^L{0) + ... (6.16) 

We thus expand all fields in the Lagrangian around their values at t = 0, making use of the equations 
of motion (6.15). The expansions up to the order we need are given by 

i^(o) = + 1(0; - yyR^.,,rx' 



e 2^ 

v;"(o) = f?" 

r{Q) = x''-{x-yyooebx' 



(6.17) 



■l,arr._im^ l {x-yy{x-yy ^a^c.b 

V (Uj - [OiUj - ^LOj (,J X 

Inserting these expansions into Sp in (6.14) yields the expression (6.13). 
Again, it is easy to check that the composition rule holds 

(ar,^|exp (-|^) Iv^x) = J d-zr^d-^.M^e-^Hx,fj\exp (^~h) exp {-^H^ \y,x) 

(6.18) 
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Imposing this rule as a consistency condition fixes the term proportional to the Ricci curvature 
times {y — xy{y — xy in the expansion of (6.11), whereas it leaves arbitrary all curvature terms 
which have a factor e in front. 

Also for this N = 2 supersymmetric case one can derive the resulting transition element from 
a path integral. As we explained in section 5, the Hamiltonian to be used in the path integral is 
the supersymmetric one, leading to the following configuration space path integral 



tation relation {'ip"' , ij/-'} = h6"''\ The last term in (6.19) is the extra term one always needs when 
writing down path integrals for coherent states as is discussed in [5] . Wc indeed find that it must be 
included to obtain the correct propagator. The actual computation of this path integral is worked 
out in the appendix. 

We will now show that the final result for the transition amplitude can again be written as the 
product of three factors: a term containing only the scalar curvature which is related to the trace 
anomaly, the exponent of the classical action, and the square root of in this case the supersymmetric 
generalisation of the Van Vleck determinant. The latter is defined by[7,21] 




where we have now rescaled the fermions by a factor ^/h, in order to obtain the usual anticommu- 




'AB 




(6.20) 



where = {x'^^rf) and = (?/■', x'')) for i^B we substitute the expressions (6.12) 

and (6.13), with the fermions rescaled by a factor Vh. To evaluate Dg write 




(6.21) 



We find, expanding in normal coordinates around x to simplify the expressions, 




(6.22) 



We do not need terms of order e in B and C, since Dg = dct^det~^(D — CA~^B^ and A^^ is 
already of order e. Writing A^^j = ^gj^f^l^d'^j + ea^j) and = ('^ab+'^aft)' write the expansion 

of the super Van Vleck determinant as 
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Multiplying by g ^/^{x)g ^^^{y) to transform D^J"^ into a bi-scalar, we obtain 

bf = g-'l\x)D''J'g-y\y) = {e)--l' 



So indeed we can write 



(x,r?|exp (-'-h) \y,x) = {2T:h)-^'^Dy\^^ ^-^(5^ + 5^)^ 



1 - ^UR 



(6.24) 



(6.25) 



similarly to (4.17). All terms involving the Ricci curvature in (6.11) are thus completely accounted 
for in the super Van Vleck determinant! This would not have been the case if we would have used 
the ordinary determinant. Note that it is Sp and not Sp which appears in the exponent, for reasons 
explained in (3.18). As we already mentioned, also in the path integral we shall need Sp. 

7. The N=l supersymmetric sigma model. 

For the case of A'^=(0, 1) supersymmetry we start again with the supersymmetric sigma model 
as defined in (5.4), but now truncate the theory by putting in the action ■0" = = 
This then requires + 62 = 0, see (5.2). Canonical quantization then yields {'0/)V'/} = 
Pi = 9ij^ + ^'^iab'^i'^l^r '^^^ truncation of would yield again (5.15) but now with ip'^i/j'^ipjipj 

instead of (V'aV'a) {i'p'^'fj ■ Using the cyclic identity of the Riemann tensor, one obtains 

-"Ar=2,truncated ~ oy 9 T^j9 Ifi"'-" 

However, this is noi the Hamiltonian we are interested in for two reasons: {i) we want a Hamiltonian 
which corresponds to the regulator of a spin-| field in quantum field theory. The latter is given 
by —^h'^g^/^lplp g^^l"^ which differs from the truncated Hamiltonian because it has a term —^fi^R 
instead of —j^h^R, and (ii), we want the A^=l supersymmetric quantum Hamiltonian, given by 
the square of the A^=l supersymmetry charge. Truncating the A^=2 supersymmetry charge in 
(5.10) to the A^=l case, one obtains the operator Q^"^ = g^^^el^{(l))ipjTrj^g~^^^. Since this operator 
corresponds to the Dirac operator g^^'^lp g~^^^ whose square contains a term \R we clearly obtain 
the Hamiltonian with —^h^R. 

The truncated Hamiltonian is thus no longer supersymmetric. The reason is that the field 
ipi — il^2 can consistently be set to zero at the classical level, as its classical supersymmetry variation 
vanishes, but at the quantum level the anticommutator of this field with itself is nonvanishing. So 
instead of the Hamiltonian in (7.1) we should use the Hamiltonian 

H = lg~'/\g'/'g'^7r^g-'/'-h'R 

in lb 
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We can now repeat the computation of the transition element of section 6. Note that in contrast 
with (6.2) we have only one species of fermions present, namely ipf. In order to compute the 
transition element (and subsequently if needed its trace), we need two operators V'" and V'" as 
in (3.16). Hence we enlarge our Hilbert space by introducing a second set of free fermions, 1^]^, 
satisfying {V'/jjV'//} = ^""^ ^^nd {^f,V'//} = 0- Wc then again define ip"" = :^(V'/ + ^"V^//) a-i^d 
i/;" = "^(V'/ ~ ^V'f/)- We rewrite the operators tpj appearing in the Hamiltonian (7.2) in terms 
of the operators V'" and '0" as ip] = -^ii^"" + V'")) and proceed analogously to the computation 
in section 6 to obtain the transition element. Compared to this section the only difference is that 
we replace the operators ^^jabCV'a = '^iahVi^^ of section 6 by icjj^j,^?^/ = I'^iahiV + + 
ip^), and Rabcd'^a'^i'^p'^p by R. The net result is that (i) aU fj"^'' terms in (6.7) and (6.8) are 
replaced by {{f]'' + r)(??^ + ^^), (n) in (6.8) the coefficient (f?"^^ry'=^'^ + 5^'=^«^<^) of the lo^ term is 
replaced by ^ [{fj"- + ^"■){ff' + ^^)(fj'' + D(^'^ + ^'^) + 2(5^''5"'^] (we have used that the expressions 
are antisymmetric under a <^ b or c <^ d, and symmetric under ab ^ cd) , and (iii) all contributions 
from the Rip^ term in the Hamiltonian disappear and are replaced by heUR. One finds 



(x,7?|exp \y,x) = (27re7l)-"/2exp (^-1(5^ + 5^) j 



1 + ^ehR{x) 



(7.3) 



where 



^ 4 (7.4) 



and is defined in (6.12). 

The various terms in (7.3) and (7.4) can be understood as follows. The fjx term comes from 
the inner product of coherent states in (6.10). The rest of the terms in (7.4) can be obtained 
from (6.13) by making the replacement fj"'x^ ~^ \{v"' + x"'){v'' + X^)- Since this substitution 
yields a result antisymmetric in a <-> 6, the (x — y^uj'^fjx term can not have a counterpart in 
(7.4). However, there is a (x — y^uj'^ term in (7.3), and this term is the counterpart of the 
{x — y^u'^fix in (6.13) as one may trace by recalling that (f/",^^ry'^^^ + was replaced by 

JE [(??" + D(^* + + + S."^) + 25'''=^'^] . The first part of this term is the square of the 

{x — y)uj{f] + x)^ term in Sp, whereas the second part explains the presence of the (x — y^uj^ term 
in (7.3). Finally, the —j^ehR from the bosonic case in (4.10) combines with the ^ehR from (7.2) 
into the term ^efiR. 

We will now show that the expression for the propagator can again be written as the product 
of the super Van Vleck determinant, the exponent of the classical action, and a term involving 
the scalar curvature which, as shown at the end of this section, determines the trace anomaly of a 
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spin-^ field. Defining the super Van Vleck determinant as in (6.20) and (6.21), we find from (7.4) 



^i6 = -^^m6W(^"+X") 

Dab = Kb + 7^{y- xyu^i^bix) + j{y- x)\y - xydiUj^^{x) 



(7.5) 



-Sconf = jdt {^^Qijx'x^ + '^s^ri^\ + '-5^r2i^i + "-x%^ri^\ ] - s^.n^^im (7.8) 



2 

which yields, using again (6.23), 

dT = (e)""/' [i - ^Riji^)iy - ^Tiy - + ^iy- ^)\y - xyuj{x)u^,'^{x)\ (7.6) 

So we can indeed write 

{x,n\ exp [~h) \y,x) = {2T:n)-^I^Df exp (^-\{Sb + )) [l + ^a'^R] (7-7) 

similarly to (4.17) and (6.25). 

Again, the same result for the propagator can also be obtained from a configuration space path 
integral with a suitable action. As exlained in the appendix, also in this case we have to introduce 
an extra set of free fermions, ip2^ in order to make the path integral well-defined. The action one 
should use in the configuration space path integral is (we again rescaled the fermions by a factor 

Vh) 

\5abri^\ + \5abr2i^l + \ 

where compared to the naive action corresponding to the Hamiltonian in (7.2) the ^fi'^R term is 
cancelled because of Weyl ordering of the bosonic term (in this A?^=l case there is no Ril)^ term 
which would yield another —\R term in the action, upsetting the cancellation of R terms in (7.8)). 
We have introduced an extra set of fermions that do not couple to any of the other fields, this way 
making certain that we do not alter the dynamics. In order to keep local Lorentz invariance we 
require that the fermions ip2 are inert under local Lorentz transformations. 

We can now easily compute the trace anomaly for a spin-^ field by taking the trace of the 
transition element (7.3) and singling out the e-independent part. Since by introducing the free 
fermions %l)2 {a = 1 . . .n), we have (for even n) 2"/^ states in the sector and 2"/^ states in the 
■02 sector (combining the n into half as many pairs of creation and absorbtion operators, and 
similarly for tp2)- Hence, we must divide the trace over tpi and V'2 by a factor 2"/^, since we really 
should only take the trace in the ijji sector. The trace anomaly for a spin-i field in n dimensions 
is therefore given by (cf. [4], equation (2.9)) 

Xr^'' =-^^mJd^Xd''fie~^''{x,fi\ew{-^H)\x,x) ; a=l...n (7.9) 

For n = 2 the trace anomaly becomes ^R, which is indeed the result for a Dirac fermion; for the 
anomaly of a Majorana fermion we have to divide this expression by two. 
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8. Conclusions. 



We have considered the evaluation of the propagator in quantum mechanics for a point particle 
in curved space, and for a point particle with its fermionic extension in curved space. In the latter 
case we considered both N=l and N=2 supersymmetric quantum mechanics. We first used the 
Hamiltonian approach and obtained the propagators through order At = e. In this approach there 
are no amibiguities at all, but the calculations are somewhat tedious. Afterwards we considered the 
problem of finding an action for a path integral, together with a presciription for evaluating this 
path integral, which reproduces the propagator. This problem has been studied at great length in 
the literature, but in our opinion no complete solution has yet been found. We now summarize the 
contribution of this paper to these problems. 

In the Hamiltonian approach, the propagator is defined by {x\ exp (^~j-H^ \y) where we assume 
H to be an operator with a given, a priori fixed, ordering of the operators and pj but without 
an explicit time dependence. We evaluated this matrix element, following Feynman, by inserting 
a complete set of states fd'^p\p){p\, expanding the exponent, and moving in each term all x to 
the left and all p to the right, where we replace them by their eigenvalues x and p, and finally 
reexponentiating the result. It seems widely believed that for small e it is sufficient to only retain the 
term linear in H in the expansion, to define {x\H\p) = h{x,p){x\p), and to replace {x\ (l — \p) 
by exp (^—^h{x,p)) {x\p). If H{x,p) is of the form T{p) + V{x), this procedure is correct, and 
proven by the Trotter formula (see e.g. [15]). However, for more general systems, such as a particle 
in curved space with Hamiltonian 

where J d"x 'y^/g{x)\x){x\ = 1 , or a point particle in flat space coupled to clectromagnetism with 
Hamiltonian 

H = I (p^ - U.ix)) (p^-iA\x)) (8.2) 

where / cPx \x){x\ = 1 , this linear approximation is incorrect. In both cases, the commutator of p^ 
with a function f{x) is given by | (^§^^ (^) (we never need to use the fact that for hermiticity p^ 
must be represented in curved space by jg~^^'*'{x)dig^^^{x) if one uses an inner product defined by 
/ (Px ^/ g{x)\x) (x| = 1 , see section 2). The reason that terms linear in e in the propagator also come 
from terms nonlinear in e in the expansion of exp (^—^H^ is that the ^-integration uses a Gaussian 
integrand exp(— e^?^), so that under the integral p is of order e~^/^. Hence each commutator 
[x,p] which removes an operator p, contributes a factor e^/^ to the final answer. Consequently, to 
obtain the terms of order e™/2 [^i the propagator for (8.1) it is sufficient and necessary to retain all 
terms with at most m commutators in the evaluation of the matrix clement (a;| ^ ^ (~t)'^ H^\p). 
In particular, terms proportional to e (which are needed to check that the propagator satisfies the 
Schrodinger equation to lowest order in e) come from terms containing none, one or two commuators 
in the evaluation of {x\H^\p) for any k\ Forgetting these commutators leads to an incorrect result. 

We have been able to resum this infinite series both for the bosonic and for the supersymmetric 
cases. It should be stressed that all terms in {x\H^\p) are well defined, and no amibiguities exist 
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in the evaluation of the propagator in the Hamiltonian approach. Only the commutation relations 
between p and x were needed, and no regularization of products of operators at the same time was 
used. Our final answer for the bosonic point particle in curved space agrees with the literature, 
proving that it was indeed necessary to keep terms beyond the linear approximation. The result 
through order e is a product of (i) the classical action for a solution of the classical equation of 
motion for x^(t), describing a particle moving from y to a; in time e, (ii) the square root of the Van 
Vleck determinant, which is anyhow needed to make the propagator a general coordinate scalar 
both in X and in y (a bi-scalar [20]), and (iii) the trace anomaly. Namely 

{x\ exp \y) = {2nh)--/'eKp (^-^S,,ix,y;e)^ e'^/^ (^1 - ^R,^ix){x - yy{x - y)^) 

exp ^-^nR{x)e^ 

(8.3) 

We would like to remark that the extension of this result to include coupling to an electro- 
magnetic field can easily be deduced from its supersymmetric extensions we considered in section 6 
and 7. One simply identifies ^A^{x) with ^^iab'^a'^Pa section 6, or with ^^iabi'i'^i section 7, 
and retains only the terms which are not coming from either anticommuting the fermions or from 
the inner product of the fermionic coherent states. This yields again the expression (8.3), where 
now is the expansion through order e of the classical action of a point particle in curved space 
coupled to electromagnetism. 

Although in (8.3) the term R{x)e can be rewritten as an integral of a local funtional of x{t), 
R{x)e = J dt R{x^i{t)) to linear order in e, the term with R^j{x){x — y)*(a; — yy can not be ob- 
tained by expanding the integral of a local functional. For this reason, one cannot simply string a 
set of matrix elements {xj\ exp (^—j^H^ l^^i+i) together and, by integrating over the intermediate 
points Xi . . .,TjY_i, obtain a path integral. The path integrals we arc interested in have, by def- 
inition, a local action, and the problem one is faced with is to find a local functional S^^^^^[x{t)\ 
which produces the propagator {x \ exp i^—j^H^ \y) when inserted into a path integral of the form 
jDx{t) exp {—f^S^^^f[x{t)]y The paths all run from x{t=—e) = y to x(t=0) = x. Of course, one 
should also specify how to evaluate the path integral, and this we do by "mode-expansion" [1,3,4], 
according to which x{t) is expanded in terms of a complete set: x{t) = J2 "n^^n(*)- The obvious and 
most convenient choice for the complete set of (f'^i't) are the eigenfunctions of the kinetic operator, 
because in that case the kinetic terms become a diagonal quadratic expression in terms of the a„. 
The problem is therefore to find a local Lagrangian i^conf [^(*)]) which can be split into a kinetic 
part -^conf and the rest, denoted by -^*Q*f, which will provide the vertices. As shown in [4] the 
correct trace anomaly, in d = 2 due to the terms proportional to e, and in d = 4 due to the terms 
proportional to e^, is obtained by using the covariant action 



1 ■ ■ 

provided one uses normal coordinates. In this paper for supersymmetric systems, and in [3] for 
a scalar particle, it was found that the complete propagator through order e is produced by the 
covariant Lagrangians we have discussed, again provided one uses normal coordinates. 
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It is not sufficient only to know fiow to calculate in normal coordinates, because some Hamil- 
tonians one encounters are not generally covaxiant. We have in mind the consistent anomalies in 
quantum field theories which maintain local scale (Wcyl) invariance and consequently break gen- 
eral coordinate (Einstein) invariance. For these systems, the regulator can be constructed using 
the algorithm of reference [12], and, as expected, the corresponding Hamiltonians are not generally 
covariant. Thus, one must know the action ^^.^^f and the rules to evaluate path integrals in arbi- 
trary coordinates. Of course, for generally covariant Hamiltonians, the propagator is a coordinate 
scalar, hence one may evaluate it in normal coordinates to obtain the correct answer. 

For general coordinates, the form of L^^^^ is unknown. In the literature, some attempts have 
been made to deduce H^^^f (related to L^^^^ by a Legendre transformation, both are functions, not 
operators) from by a functional differential equation. In particular, DeWitt has pursued this 
problem in the latest version of his book on supermanifolds[7], but he runs into the well-known 
problem how to define products of operators at the same point (i.e., at the same time). He proposes 
to use "time-ordering", according to which one first defines all operators at different times, then 
considers all permutations of these operators, and finally takes the limit that all times become 
equal. Following this procedure for the bosonic case, one obtains not only the term mentioned 
above, but also an extra noncovariant term 

^conf = Ig'^PiPj - l^'R - Ih'T^jkr'ug'^ (8.5) 

Concluding that time-ordering breaks general coordinate covariance, he chooses a pseudo-Euclidean 
frame in which r*^-^(x(t=0)) = and computes (for the supersymmetric case) some two loop 
diagrams and finds that they are in agreement with the result for the propagator (which he obtained 
long ago by heat kernel methods). 

We believe that this is not the whole story. First of all, the presence of a F^ term at all t will 
show up in calculations using normal coordinates at order (where 3-loop diagrams contribute). 
Although these contributions do look covariant in normal coordinates, one obtains an incorrect 
result for the d = 4 trace anomaly. Hence the ^F^ term without further modifications should not 
be present in L^^f when one uses normal coordinates. Of course, once the correct action S^.^^^ 
is known in normal coordinates, one can find the corresponding action in general coordinates by 
making a change of integration variables. The Jacobian will contribute new vertices, and although 
we have not worked out the path integrals for general metrics, we have a remark. The background 
decompositon x*(t) = XQ(t) + a^qu(i) seems too linear to us for such a nonabelian theory as gravity. 
Let us recall that in superspacc Yang-Mills theory, the Yang-Mills field is contained in exp{V), 
and a background field approach replaces this by exp(Vback) ^^P(^u)- This is clearly an infinitely 
nonlinear split. Analogously, we expect also in gravity a nonlinear split of the form x'^{t) = Xg(t) -|- 
^qu(0 + (iCqu(^))^ terms -|- . . .. We do not claim to have exhaustively studied this problem , and 
intend to return to it in the future. However, we have made some introductory calculations which 
we now describe. 

In our explicit computation of the path integral to order e (2 loops), we took Lq = ^g^j{x)x''x^ 
where the x in g^j {x) is the endpoint of the interval. Decomposing the action into the sum of this 
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kinetic term and vertices we computed loops. We then proceeded to evaluate all diagrams up to 2 
loops in an arbitrary coordinate frame 



^ = 11'^^ Idijixoir) + Xqjr)) [(4(r) + ir^Jr)) + i'qu(^)) + 6'(r)c'(r) + a\T)a^ {r) 

(8.6) 

where t = t/e. Here a, b, c are anticommuting ghosts due to exponentiating the measure. They 
were first introduced by Bastianelli in [3]. The present form was obtained in [4]. Due to these 
ghosts all loop contributions become completely finite. We contract all quantum fields, using (see 
(A.7)) 

< <u(f^Xu(^) > = -eng'^A{a,T) 

< b'{a)c'{T) > = -2ehg'^dlA{a,T) (8.7) 

< a\a)a^{T) > = ehg'^ dlA{a,T) 

where 

A/^^^ sin(n7ra) sin(n7rT) 

A{a, r) = -2 2^ (8.8) 

n=l 

The great advantage of the formulation in (8.6) is that it allows a regularization which treats x and 
the ghosts in a uniform manner, namely by mode cut-off: in all graphs one first uses A(a, r) in (8.8) 
where the sum runs up to N, and only at the end takes N to infinity. Of course, Xo(i) = x+{y—x)t/e 
is only a solution for Sq and not a solution of the full equation of motion, hence there are terms 
linear in Xq^{t) in the action. The path integral expectation value of a single Xq^{t) vanishes (since 
Lq is symmetric under Xq^{t) ^ ~x^^{t)), but contractions of the form x^^{s)x^^{t) contribute. 

We found that the total set of 2 loop diagrams in the bosonic case precisely reproduces all 
terms in the propagator through order e, except for one term (which vanishes in normal coordinates). 
Namely, if one uses the covariant action (omitting the term), the path integral result is equal 
to the propagator plus the following remainder 

-^ehg'^g'^^g-^- {d,g,J {d^g^^) + le^^^^/'^-" {d,g,J (digj,) = -Lehg^.V^.^r^ / (8.9) 

If one adds the term, part of this remainder is canceled, but we are still left with another extra 
term 

{d,g,J {digj^) (8.10) 

We believe that a proper treatment of the Jacobian for the transformation from normal to general 
coordinates will cancel (8.9). Our conclusion for the bosonic case is that in the Hamiltonian sector 
we have a complete understanding, but in the path integral sector, we believe we know the correct 
'S'conf (namely the covariant action with the term but without terms) but we have no proof 
of this rule. 

We now turn to the fermionic sector where we have not only found similar results, but also 
new results for fermionic phase space path integrals. In the N=2 supersymmetric case, there 
are also fermionic variables tpai't) present, with a = 1,2 and a = l...n. The Hamiltonian for 
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this system is not obtained from the regulator of a corresponding quantum field theory, but we 
fixed its ordering by requiring the supersymmetry relation 25^pH = {Qa^Q/s} where the operator 
orderings in the supersymmetry charge operator in (5.10) were fixed by requiring hermiticity. 
We then computed the propagator between states \y,S,) and {x,f]\. The states |,^) and (f^| are 
fermionic coherent states, and using general theory of coherent states, suitably adapted to the 
fermionic case, we straightforwardly obtained the propagator to order e. Of course, we must now 
also take into account the commutators {V'^", ^/'^j = S"-^. The final result in (6.11) was then written 
as exp {—-^{Sb + Sp) — (^ab^^^M^)) times a factor containing only terms with curvatures. The 
extra term —Sai,fj""tpci{^) (where V'ci(i) is a solution of the field equations with boundary value 
V'ci(i=— e) = X^) is well-known from the theory of path integrals for coherent states, and can be 
found in textbooks, see section 3. The curvature terms should describe the Van Vleck determinant 
and a term with R which would have been the trace anomaly of a quantum field theory if H could 
have been identified with its regulator. However, in supersymmetric theories it is more natural 
to use superdeterminants, and hence we considered the double super derivative of the action. The 
super Van Vleck determinant absorbed then all terms with Ricci curvatures, just as in the bosonic 
case. We were left with a term involving the scalar curvature, which is exactly the same as in the 
bosonic case. (Using the ordinary Van Vleck determinant, one is left with a term containing the 
Ricci curvature). Thus, the propagator factorizes exactly as in the bosonic case, but now with a 
super Van Vleck determinant. 

To reproduce these operator results for the N=2 supersymmetric case from a path integral, 
we took as action S^^^^f the classical N=2 action in (6.19) with the boundary term, including the 
Rip'^ term but without any further term proportional to R. The reason we did not include an R 
term is that, applying the time-ordering prescription to the N=2 Hamiltonian, we have found that 
the term — -^i? of the bosonic action is precisely canceled by a similar term obtained by time- 
ordering the Rip^ term! (In the N=l case a similar cancellation of |i? terms occurs for different 
reasons). Time-ordering however, also produces the term of the bosonic sector discussed above, 
together with an term if one treats p, x, -ip and as the independent variables. (However, with 
— u^gjj{x)tl) tp°, X, i/j and -0 as independent variables, no lo term is produced). 

With the fermionic propagator given by the mode expansion (and with the propagator for 
the bosonic and ghost fields mentioned before) we then may compute all 2-loop graphs using (8.6) 
and (8.7) together with (A. 15) and (A. 16). The result of these calculations can be summarized 
as follows: to order e, the path integral produces the propagator plus two extra terms. The first 
one comes from the bosonic sector and is equal to (8.9), the second one is due to contracting two 
Jc/Ta;q^a;jg^(a;)V5q„V'qu vertices. Writing the three propagators in a mode expansion, and truncat- 
ing the bosonic and fermionic propagators at the same n = N, we have found (by a numerical 
evaluation) a finite but nonzero result 

-^e%,,r^^rV - ^ehg^^cojco^,^ (8.11) 

In normal coordinates uj^^f,{x) vanishes, and hence in normal coordinates we obtain complete agree- 
ment between the operator approach and path integrals. In a future extension to general coordi- 
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nates, the contribution in (8.11) should be taken into account. 

Finally, we considered the supers3mimetric N=l case. This case is more complicated than the 
N=2 case, because there are only real fermions present, whereas one needs complex fermions, both 
in the Hamiltonian approach (for operators and ^^^) and in the path integrals (to define separate 
boundary conditions on the left and the right). We resolved this problem by adding a set of free 
fermions in both cases; this should not alter the model. To preserve local Lorentz invariance, we 
required that these free fermions are inert under these transformations. 

In their fundamental paper on chiral and gravitational anomalies [1], Alvarez-Gaume and Wit- 
ten did not need to introduce these extra fermions because they only considered traces of the 
propagator, so they did not need to specify separate boundary conditions on the left and the right. 
Moreover, the Jacobian for chiral transformations of the quantum fields led in the corresponding 
quantum mechanical model to periodic boundary conditions for the N=l case, and half-periodic 
and half-antiperiodic or completely periodic for the N=2 case. In all these cases, there is a zero 
mode ipQ in the expansion of ip{t) in terms of cigcnfuncions of the operator ^ (corresponding to our 
Sq). The Grassmann integration then requires a minimal number of eRip'^ipQ vertices to saturate 
the integration over these zero modes, and this already produces enough powers of e to cancel the 
g-n/2 -j^ ^]^g measure of the path integral. Consequently, for them the harmonic approximation 
is sufficient. For us, however, this is not sufficient. For example, for the trace anomaly we need 
completely antiperiodic boundary conditions for the fermions, so that no e-producing zero modes 
ipQ are present. In addition, wc considered the full propagator, and not only its traces. 

The actual computations in the A^=l case became very similar to those of the N=2 case, 
once the free fermions were added. The Hamiltonian which corresponds to the general coordinate 
invariant regulator of a spin-^ field, was the sum of the Laplace-Beltrami operator for fermions 
plus a term —^h^R (twice as much as truncation of the N=2 Hamiltonian would give). This result 
was consistent with the requirement that H be equal to the square of a hermitian supersymmetry 
charge. Again we found that the propagator factorized into the classical action (with boundary 
term), the super Van Vleck determinant, and a scalar curvature term which yields now the trace 
anomaly for the quantum field theory of a (Dirac or Majorana) fermion minimally coupled to 
gravity. 

Then we considered the path integral approach to the propagator of the N=l case. In this 
case we took as action S'^.Q^f the action in (7.8), which is the classical supersymmetric action (with 
boundary term). No extra R term was present because wc showed that the R term due to time- 
ordering in the bosonic sector cancels the explicit R present in the Hamiltonian. Hence, for different 
reasons, both in the N=l and N=2 cases, the action S^^f equals the classical supersymmetric 
action. This result one might have taken anyhow on the basis of symmetry arguments, but it 
is gratifying that if follows rigorously from our Hamiltonian approach. We computed all graphs 
through two loops which contribute through order e. In normal coordinates we found again complete 
agreement with the operator approach, but in general coordinates we found again an extra term 
which a future extension of the path integral approach using arbitrary coordinates should explain. 

We conclude this paper with a general solution for the path integral in terms of an action 
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where g = det^j^. Then use the metric gij{x) to construct normal coordinates. In these normal 
coordinates one has 



The path integral is evaluated by splitting off a kinetic term g^j{x{0))x^(t)x^ (t), and adding the 
ghost action in (8.6). Propagators and vertices are now defined, and the loop expansion of the 
propagator in normal coordinates is obtained. Using mode cut-off [3,4] to regulate divergent graphs, 
it is found that at any loop level the final result is completely finite. For general coordinate invariant 
'S'conf' t'^is result also holds in general coordinates. In other cases one must transform back from 
normal to general coordinates. We have checked that this procedure reproduces the propagator 
obtained from H through order e, and trace anomalies through order e^, for the point particle, 
and A^=l and N=2 supersymmetry. Summarizing: whereas the heat kernel or the Hamiltonian 
approach yield the propagator directly, we have found the corresponding action S^^^^ which exactly 
reproduces this propagator from a path integral. 
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A. Appendix: Evaluation of path integrals. 

The path integral in the bosonic case (section 4) has already been evaluated in normal coordi- 
nates in [3,4], by expanding fields in Fourier series and converting the path integral to an integral 
over the Fourier coefficients. An important detail is that we must introduce a ghost system to ex- 



ponentiate the factors ■^g{xj) at the N — 1 intermediate points Xj in the measure, hereby maJiing 
the remainder of the measure translationally invariant. This makes loop corrections finite. We will 
now first evaluate the bosonic path integral in an arbitrary coordinate frame, and afterwards we 
will compute the fermionic path integrals. 

We follow the approach of [3,4], and start with the action in configuration space (note that 
extra terms of the type ii + are already of order e, and can, to the order we are considering, 
simply be replaced by their classical values) 








(8.13) 



-At 









(A.l) 
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where we have rescaled t = t/e, and introduced a ghost system to exponentiate the factors g{x{t)) 
appearing in the measure. The boundary conditions arc given by a;*(— 1) = y*, a;'(0) = x*. Ah ghosts 
vanish at both endpoints, 1) = 6*(0) = c*(— 1) — c*(0) = a*(— 1) = a*(0) = 0, because there 
are only factors at the intermediate points. We will now make a decomposition in background 
and quantum fields, where wc choose the background fields to satisfy the boundary conditions. We 
take vanishing background fields in the ghost sector, whereas for the bosonic fields we take 

x\t) = xi{T) + x^^^{t) ■ x'o(r) = x^ + (x-y)V (A.2) 

We now evaluate the path integral by expanding the term gij{x{T)) around its value at x(t) = x, 
and write S = S'^'" + S*"*, where 



gkin = ^JdT [x'x^ + 6V' + a'a^] (A.3) 
-1 

and 



5-* = ^ /dr + \d,dig,^x''x'^ \x'x^ + b'c^ + a'a^] (A.4) 

-1 

Now all functions g, dg are evaluated at the endpoint x. Since all quantum fields vanish at r = and 
at r = — 1, we can decompose them in a mode expansion, yielding a natural regularization scheme 
for the propagators. The propagator was obtained in reference [4] by using the mode expansion, but 
one can also obtain it (more easily) from canonical methods. Expanding the Heisenberg operator 
X in a complete set of zero modes of the kinetic operator ~5'tj(-^)^ we get x*(t) = x* + g^^ (x)pjt, 
with [x^,Pj] = ihdj. Imposing the boundary conditions Xqu(t=— e) = Xqu(t=0) = as conditions on 
the bra and ket vacua, x{t=-e)\0) = and {0\x{t=0) = 0, we find (0|x* = and (x* — ep*)|0) = 0. 
Hence 

A'^{s,t) = {0\Tx\s)x^ {t)\0) 

= e{s-t){0\{x' +p's){x^ +ft)\0) + (s^t) 

= e{s-t){0\p'sx^{l + t/e)\0) + {s^t) ^^-^^ 

= e(s-t)-g'h(l+t/e) + (s^t) 
i 

The Euclidean version of this result, after we make the rescaling t = t/e, is given by 

A'^ {a, t) = -ehg'^ [a{T + l)9{a - t) + T{a + 1)9{t - a)] (A.6) 

Using the mode expansion and (A.3) we find the same propagator, but now with a regularized 
expression for A((t, r) (see [4]) 

n=l 



The simplest term comes about when we only consider the background fields Xq{t). This yields 



{y-xy{y-xy (A.8) 
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which already reproduces most terms in the classical action (compare (4.11)); note that since we are 
not expanding around the full solution to the equation of motion, we cannot expect to reproduce 
the full classical action in Sq. We will see that the remainder of the classical action is produced by 

quantum corrections. 

Next we consider contributions from terms in which we contract one pair of quantum fields. 
These come about either when we just take 5*"* by itself, or when we take (S**^*)^. The terms of 
the first type yield 



^g'^ix - yf{x - - g>\x - y)\x - yf + ^g^\x - y)\x - yY 



(A.9) 



12 

which is part of the Rij{x — y)''{x — yY term. Furthermore, the terms of the second type yield 

l^ff^r^r^K^ - yy{^ - vYi^ - y)H^ - vY (A.io) 

which we recognize as the remainder of the classical action. 

The next set of terms comes from contracting two pairs of quantum fields. Again there are 
two types of contributions, from S*"* by itself and from its square. The first terms are 

^,nd,dig,j[g''^g^^-g^''g^'] (A.ll) 

which we can identify as part of the scalar curvature term, whereas the terms coming from (5*"*)^ 
combine with the terms in (A.9) into —j^Rij{x — yy{x — yY . 

Finally, we should contract three pairs of quantum fields, now only in (5'"*) . Part of these 
terms we need to covariantize the expression in (A.ll) into ^ehR. The remaining terms are given 

by 

-^ehg^'g^'g'^'' {digkJ {djgin) + ^ehg'^g'^g^^ {d,g,J (9,5,J (A.12) 

which is equal to 

-^ehg,ir\^T^/ (A.13) 

The significance of this term is discussed in the conclusions. 

We now evaluate the fermionic (supersymmetric) path integrals. In these cases, the bosonic 
and ghost sector of the path integral can be handled in exactly the same way as before, again 
yielding the action together with the contribution ^ehR — ^Rj^j{x — yy{x — yY . We now 
consider the fermionic sector. 

Similarly to [3,4], we rescale r = t/e to make the e-depcndence more explicit and facilitate 
keeping track of the order in the expansion in e. We now first note that we can write 

S = S^ + S^:^ + Sl^i (A.14) 

Here contains all terms with only bosonic or ghost fields. The other two terms are given by 



S^ = JdTS,,r^^ (A.15) 
-1 
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and 



Qint 
'-'fer 
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-I 

-1 



dr 



1 



(A.16) 



for the N=2 supersymmetric nonlinear sigma model considered in section 6, whereas for the N=l 
model discussed in section 7 we have 
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■I 

-1 



dT 



(A.17) 



In a background field approach we decompose V'(r) = '4^q{t) + V'qu (''")> = '0o(''") + 

?^qy(r). We could have used fermionic normal coordinates and added a term proportional to x^^ as 
in [2], but the result should not depend on the split between classical and quantum fields, and our 

choice is quite simple. We choose V'o(t) and V'o(''") to be the solutions the equations of motion of the 
kinetic part of the action, (A. 15), which satisfy the boundary conditions, i.e. we take '0q(t) = f)"" 
and 'iPq{t) = x". Since wc require ^"(0) = fj'^ and ip"'{—l) = x", this implies that the quantum 
fields need to satisfy the boundary conditions V'qu(O) — V'qu(~l) — 0- The Green function 
for the action (A. 15) of the quantum fields then reads 



G%{a, t) =< Jc7X,^(r) >= -M«^,5,^(0(a - r) - 0(r - a)) + 



yf3 



(A.18) 



where = +1- Or, in terms of and ip, 

< ^ui^Xuir) >-- 



ns^^eir - a) 



< r^ni^Xnir) >= =< <(a)V;^Jr) > 



^qu^ 



(A.19) 



One can easily check that this Green function is a solution to the equations of motion of (A. 15) and 
that the boundary conditions are indeed satisfied. In particular, at cr = r we define < A{t)B{t) >= 
lim^^o 5 [M'T + S)B{t - S) + A{t - 5)B{t + 5)] , and at(7 = r = 0orcr = r = -l one gets the 
correct result using ^(0) = \. 

The propagator for the fields V'(i) and tp^t) can also easily be found from the mode expansion 
or from canonical methods. In the canonical approach -(/'(t) = ip and ^(t) = ^ with boundary 
conditions on the vacua given by (0|'i/'(i=0) = and il}{t=—e)\Q) = 0. In the mode expansion, we 
write 



cos 



((n+-)7rt/e) 



n 



2' 
1, 



(A.20) 



and the kinetic matrix J dt^V" becomes ^ \T^{n + 5)^71 yielding the propagator 

— e 

< T{[)''{s)i)\t) > 



n=0 ^ 



Y sin((n + l:)TTs/e) cos((n + ^)7ri/e) 



(A.21) 
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This is the ^-function < T-ijj"' {s)'tp^ {t) >= 6"'^6{s — t) as follows from a Fourier analysis. For our 
purposes we prefer the expression (A. 21) for the propagators in terms of modes, since this yields a 
convenient regularization. 

We can now straightforwardly compute the transition element by expanding exp {—%Sl^^^ and 
contracting the quantum fields, using the above Green functions for the fermions. For the bosonic 
and ghost fields, the Green functions have been computed in [4]. The propagator for the bosonic 
fields is given in (A. 6). 

We will first consider the N=2 case, with the interaction given in (A. 16). When we expand 
exp {—%Sl2r) we will for the first term in this expansion only need the contraction at equal time 
< '0"(t)'0^(t) >= —^^(5"^. Since this is symmetric in a and 6, the first term in (A.16) will yield 
no contribution, whereas the second term contributes ^ei^^^ry^x'' — ^ehR. Next consider the term 
I {h^icr)"^ ■ Terms involving the term are of higher order in e (for example, the cross term is of 
order (x — y)e ~ e^*^^), so we only need to find the contribution from the square of the xujip'^ term. 
When we contract only the four fermionic fields, using (A. 19), one finds zero. The contraction of 
the two X fields and two fermionic fields (using (A. 19) and (A. 6)) vanishes because the fermionic 
contractions yield 6{a — T) + 9{t — a) = 1, and the integral of the bosonic part vanishes. Finally, we 
can contract four fermionic and two bosonic fields, which yields a contribution —^efig^^Uj^J'ojj^^. 
In order to obtain this result, we used the regularized version of the propagator from the mode 
expansion, evaluated the integrals using the same number of bosonic and fermionic modes, and 
only afterwards let the number of modes go to infinity. Contractions involving the other bosonic 
fields are again of higher order and need not be considered. 

We should be careful since we are in fact expanding the fermionic fields V'"(r) and ijj'^{T) 
around the constant fields 7)°" and respectively, which are not solutions to the equations of 
motion (6.15) of the full action. Since < V'qul''") >=< V'qul''") *-* due to the fact that the action 
is even in the number of Grassmann fields, the term from Sl^^ which is linear in quantum fields 
vanishes. However, we do have to take into account the contribution when we contract only one 
pair of fermions, each field from a different factor Sl^^ in ^ {\Sl^r)^ ■ This yields the contribution 
|(x* — y^){x^ — y^)uj^g^'^u}j^f)^-)^ . When we were expanding around the solutions to the equations 
of motion, this term was part of the classical action (6.13). Since we are now expanding around the 
constant fields f)"" and x", this term is not part of the expansion of the action but instead appears 
as an additional quantum correction. 

Taking the contributions from the bosonic and ghost sector, together with the above computed 
contributions from the fermionic sector, we arrive at (6.11), plus a remainder — j^ehg^^ cOj^J'LOjjj"' — 
^efigj^lT^j^T^ . Similarly to the bosonic case, this term vanishes in normal coordinates, so that 
we obtain complete agreement between the Hamiltonian and path integral approach. In general 
coordinates, we again expect that a proper treatment of the Jacobian generated by a change of 
variables from normal to general coordinates will lead to cancellation of this term. 

Finally we consider the A'^=l case. Our aim is to reproduce the result in (7.3) from a path 
integral with (7.8) as action, and with boundary conditions V^"(0) = ry" and ip"'{—l) = x". As 
explained in the beginning of this appendix we should obtain from the fermionic sector the terms 
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in Sp in (7.4) but also the last term in (7.3). In this case it is easier to work with the ipi propagator 
given in (A. 18) then with the -ip, ^ propagators, since the fermionic sector of the interactions, given 

in (A. 17), only depends on ip^. We again decompose ^pf into a sum ijjf^Q + V'l.qu with ■0"o = 
■^j^iv"' + x") ^-iid V'l qu appearing in the propagators. The first term in Sp is of course due to the 
last term in (7.8). The next set of contributions comes from the expectation value of —y^Sl^^. The 
classical part is given by 



-1 



(A.22) 



where Xq{t) is a solution of the bosonic part of the action. Expanding Xq(t) = x* + (x — y)*r + 

0{x — y)^, one recovers the last two terms in (7.4). The equal time contraction of V'lV'i i^i (A. 17) 

vanishes by itself (it anyhow is proportional to d""^). Also the contribution from the contraction in 



x'Wjjjf, in (A. 17) vanishes, as it is proportional to J dr < x*(t)x-^ (t) > which is zero according to 

-1 

(A. 6). We are left with the last term in (7.3), which should come from the term i . The 

classical part of this term is of course contained in the expansion of exp (^—^Sp^ . The contraction 
of x*((t) with x^{t) is of order e, and would contribute, but its integral over a and r vanishes. One 
x*(o") with x-' (r) contraction is also of order e, but would leave a factor Xq which is of order e^^^, 
so this term is of higher order. The contraction of all four fermions is nonvanishing, and indeed 
reproduces the term ^(x* — y^){x^ — y^)uiiJ'u)ji^°' , in agreement with (7.3). Finally, we can contract 
four fermionic fields and two bosonic fields x*. This yields — -^ehg^^ u^J'ujjij"- , where, as in the N=2 
case, we take the same number of bosonic and fermionic modes to regularize the integrals. Adding 
all contributions we find the propagator for the A^=l case as given in (7.3), plus a remainder 
which now equals —■^ehg^^uJiJ'ujji,"' — ^eTjg^.^r'^j-^T'j*. Again, for general coordinates we expect 
this term to cancel against a Jacobian; in normal coordinates we find complete agreement with the 
Hamiltonian result. 
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